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INVESTIGATION OF STEADY AND FLUCTUATING PRESSURES 
ASSOCIATED WITH THE TRANSONIC BUFFETING AND 
WING ROCK OF A ONE-SEVENTH SCALE MODEL OF 
THE F-5A AIRCRAFT 


!-• 

a 


By Chintsun Hwang and W.S. Pi 
Northrop Corporation, Aircraft Division 
Hawthorne, California 


SUMMARY 

The report describes a wind tunnel test of a one-seventh scale F-5A 
model and the evaluation of pressure, force, and dynamic response measure- 
ments during buffet and wing rock. The detail instrumentation, the specially 
designed support system which allowed the model to oscillate in roll to sim- 
ulate wing rock, are also described. The investigation was conducted under 
Contract NAS2-8734, sponsored by NASA Ames Research Center. 
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The dynamic pressure and response data were recorded on frequency- 
modulated (FM) tapes which were used after the test for real-time and spectral 
analyses. The static and rms fluctuating pressure data were reduced to non- 
dimensional coefficients which were applied to buffet intensity and wing rock 
onset studies. The dynamic pressure data obtained during the wind tunnel tests 
were compared and correlated with those previously recorded in a flight test 
program. Parametric studies were performed using the acquired test data 
(by varying either the model configuration or the test condition) to investi- 
gate the effects of control surface settings, Mach number, sideslip angle, 
etc. on buffet. 

Based on the pressure (steady state and fluctuating) and response data 
acquired when the model was near stationary and when it was excited by flow to 
oscillate in roll, a limit cycle mechanism was identified which supplied energy 
to the aircraft model and caused the wing rock type oscillations. The major 
origin of the fluctuating pressures which contributed to the limit cycle was 
traced to the wing surface leading edge stall and the subsequent lift recovery. 
For typical wing rock oscillations, the energy balance between the pressure 
work input and the energy consumed by the model aerodynamic and mechanical 
damping was formulated and numerical data presented. 
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Section 1 


INTRODUCTION 

The work described in this report deals with the wind tunnel test of 
a one-seventh scale model of the Northrop F-5A aircraft at transonic speed 
to investigate the steady and fluctuating pressures associated with buffet 
and wing rock. The tests were performed at NASA Ames Research Center 
Eleven-Foot Transonic Wind Tunnel. 

Under sponsorship of the Ames Research Center (ARC), Northrop previously 
completed Contract NAS2-6475, entitled "Investigation of Northrop F-5A Wing 
Buffet Intensity in Transonic Flight," That project was a flight test pro- 
gram to acquire and analyze F— 5A aircraft wing dynamic buffet load and 
response data in the transonic region. Based on the flight test data 
evaluation, the pressure distribution on the aircraft wing during buffet was 
used to define the initiation of the shock front(s) and the increase in 
dynamic pressure intensity due to flow separation, as well as the expansion 
of the separation region as the angle-of— attack increased. Within the tran- 
sonic .region, the effects of the Mach number with a constant-q load, as well 
as the effects of various flap settings on the pressure distribution and air- 
craft responses, were delineated and correlated analytically. The flight test 
results are described in Reference 1. In order to compare and correlate 
the dynamic pressure data acquired during the wind tunnel test with those of 
the flight test, taking into consideration the scaling relations, the 
scale model tests were conducted using test conditions similar to the flight 
test conditions. 

The F-5A scale model wind tunnel buffet investigation program was ac- 
complished in two phases. In the first test phase, a nominally rigid sting 
was used. In the second test phase, a special sting was designed incorpor- 
ating a torsional spring and damper which allowed the model to oscillate in 
roll at a natural frequency simulating the Dutch roll motion during wing— rock. 
Special instrumentation for the model included 28 Kulite Model LQ2-156W-4 
pressure transducers, 28 static pressure taps, 3 accelerometers, and 4 wing 
root strain gage bridges. Tests were performed encompassing 11 model con- 
figurations defined by the external stores and control surface settings. Three 
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Mach numbers were selected for testing which were nominally identical to the Mach 
numbers used previously in the flight test program^ Other test parameters 
were the Reynolds number and the sideslip angle. Corresponding to each test 
condition, various angle-of-attack settings were used up to +16 degrees, and 
the data were recorded for all s in 2 to 4 degree increments. 

During the wind tunnel tests, pressure, force and dynamic response 
measurements were recorded on frequency-modulated (Ffi) tapes for real-time and 
spectral analyses. The force, static pressure, and rms fluctuating pressure 
data were reduced to coefficients used for buffet intensity and wing— rock 
onset studies. Substantial amounts of tests, data processing, evaluation and 
correlation work were carried out and documented. This report consists 
of the following key work items. 

1. Study of the stationer ity of the real-time dynamic pressure data. 

2. Study of the effect of Mach number, wing— tip— mounted missiles and 
tail deflections on buffet intensity. 

3. Study of the effect of the flexible support on buffet intensity, 

4. Study of the development of flow separation bubble on the scale 
model wing with sideslip. 

5. Comparison of spectral data based on the wind tunnel test and the 
flight test. 

6. Study of the effect of flap settings on buffet intensity, 

7. Evaluation of the model response data. 

8. Investigation of the wake flow effect and tail surface buffet 
intensity data. 

9. Correlation study of dynamic pressure data under sideslip 
conditions. 

10. Investigation of the Reynolds number effect based on the wind 
tunnel test data and corresponding flight test data. 

11. Presentation of low-frequency-range power spectral density plots 

of fluctuating pressures and responses when maximum wing-rock motion 
occurred. 
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Determination of the system aerodynamic damping by applying the 
Randomdec program corresponding to the model wing-rock. 

13. Investigation of the aerodynamic hysteresis effect on F-5A wing-rock 
behavior based on a limit cycle concept. 

14. Nondimensional data presentation and comparison with similar data 
obtained in other separated flow tests. 

15. Analytical formulation of scale model modal data generation and 
stability prediction allowing for aeroelastic effect. 

The program was conducted under the sponsorship of Ames Research Center, 
National Aeronautics and Space Administration. The NASA program monitor 
is Mr. Charles F. Coe. Messrs. Donald Buell and Dennis Riddle of Ames 
Research Center participated and contributed substantially in the conduction 
of the wind tunnel tests. At Northrop, other than the authors of this report, 
the key contributors included Dr. Wilford Wong and Mr. Brent Bennett. The 
flexible support device of the test model was designed by Messrs. Mel C. 
Sanders and John E. Black. 
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Mean aerodynamic chord or wing span, meters 
Damping ratio 

Nondimensional aerodynamic derivatives 
Static pressure coefficient 
RMS pressure coefficient 
Frequency, Hz 

Generalized force matrices, kilonewtons, etc. 

Flight altitude or wing translatory movement, meters 
Reduced Frequency 
Characteristic length 
Mach number 

Total mass or total mass moment-of-inertia , kilograms, etc. 
Non-dimensional frequency 

2 

Fluctuating pressure, kilonewtons/ (meter) ; p also the roll 
rate, radians/second 

2 

RMS value of the fluctuating pressure, kilonewtons/ (meter) 

2 

Static pressure, kilonewtons/ (meter) 

2 

Total pressure, kilonewtons/ (meter) 

2 

Free stream dynamic pressure, kilonewtons/ (meter) 

Yaw rate, radians/second 

Reynolds number, (meter) ^ or non-dimensional 

2 

Wing surface area, (meter) 

Freestream speed, meters/second 

Spanwise distance from the roll axis , meters 

Angle-of-attack, degrees 

Critical angle-of-attack, degrees 

Sideslip angle, degrees 
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Aileron deflection, degrees 
Trailing edge flap angle, degrees 
Horizontal tail deflection, degrees 
Leading edge flap angle, degrees 
Eigenvalue 

Maximum correlation coefficient 

Time delay, seconds 

Roll oscillation angle, degrees 

2 2 

Power spectral density of pressure (kilonewtons/meter ) /Hz 

Circular frequency, radians/second 

Natural circular frequencies, radians/second 


Matrices 
[a],[A( )] 

L^( )] 

)1 

1 f( ) j 

[ % J 
)! 

)] 

'•^fI 

[t] 

j xj, [x] 

fzj 

)«} 


Generalized AIC matrix, newtons/meter 
Damping coefficient matrices, newton — second/meter 
Fourier transforms of force matrices, newtons 
Generalized force matrices, newtons 
Modal damping coefficient matrix 
Fourier transforms of deflection matrices, meters 
Structural stiffness matrices, newtons/raeter 
Modal stiffness matrix, newtons/meter 

Stiffness matrix for the flexible support, newtons /meter, etc. 

Rigid body transformation matrix 

Modal deflection matrices 

Modal impedance matrix, newtons /me ter 

Modal amplitude matrix, meters 
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Section 2 


WIND TUNNEL TEST PROGRAM 

The buffet wind tunnel test program was carried out on a one-seventh 
scale model of the F-5A airplane at Ames Research Center Eleven— Foot 
Transonic Wind Tunnel, The objective of the test program is to obtain wind 
tunnel pressure, force, and dynamic response measurements with which to 
develop methods for predicting buffet intensity and wing-rock instability at 
transonic speed. The program was accomplished in two entries. During Phase 
I, the first entry, a standard sting was used. The second tunnel entry (Phase 
II) was carried out with a flexible support system that allowed for model 
oscillations, especially those degrees of freedom corresponding to the Dutch 
roll motion during wing- rock. 

Previous to the model tests, a buffet flight test program was conducted 
using a fully instrumented F-5A aircraft. The test results of the flight test 
are described in Reference 1. The scale model tests were conducted using 
test conditions similar to the flight test conditions so that the dynamic 
pressure data acquired during the flight test and the scale model tests may 
be compared and evaluated, taking into consideration the scaling relations. 

The F-5A is a single-seat fighter capable of carrying stores at wing and 
fuselage pylon stations. The flight test was conducted with two wing-tip 
stores (AIM-9B Missiles) with guide rails; otherwise the wing was clean. 

The scale model tests were conducted with and without the wing— tip missiles. 

A combination of deflected leading edge and trailing edge flaps as well as 
the case of completely retracted flaps were used in the test program. A 
three-view drawing of the F-5A is shown as Figure 1. 

The basic wing airfoil of the F— 5A is a NACA 65A004.8 section modi- 
fied as follows. From the 407. chord to the 66.67. chord, the thickness 
distribution corresponds to a NACA 004.8-64 airfoil. Aft of the 827. chord, 
the airfoil is composed of straight lines forming a 7® trailing edge angle. 
Between the 66.67. chord and the 827. chord, the section has been curve- 
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Airfoil Section 
Area (Reference) 

Span (clean tlpe) 
Aspect Ratio 
Taper Ratio 
Sweepback (251 Chord) 
Mean Aerodyriaolc Chord 
Dihedral Angle 
Incidence Angle 


MACA 6SA00A.8 (Modified) 
15.79 (170.00 ft2) 

7.696 o ( 25.25 ft ) 

3.75 
.20 
24 * 

2.356 m ( 7.73 ft ) 

0 
0 




fitted so that the ordinates, slopes, and curvature at the end points are 
continuous. The airfoil is cambered to a NACA (.65) 50 five-digit series 
mean line. The mean line is rotated so that the straight aft portion lies 
in the wing chord plane and the camber appears as a 17. leading edge droop. 
The airfoil is constant over the entire span. 


The horizontal tail of the F-5A is all-movable and has an unmodi- 
fied NACA 65A004 section with the following physical dimensions and con- 
figurations : 


Area, Total 

5.48 m 

(59 ft ) 

Area, Exposed 

3.07 m^ 

(33.03 ft^) 

Aspect Ratio, Exposed 

2.88 


Sweepback (257. Chord) 

25° 


Mean Aerodynamic Chord, Exposed 

1.222 m 

(3.68 ft) 

The vertical tail airfoil is a modified NACA 

65A004 section. The air- 

is composed of straight lines aft 
) of the vertical tail are: 

of the 727. 

chord. The basic dimen- 

Area, Exposed 

3.83 m^ 

(41.2 ft^) 

Aspect Ratio, Exposed 

1.22 


Taper Ratio, Exposed 

.25 


Sweepback (257. Chord) 

25° 


Mean Aerodynamic Chord, Exposed 

1.990 m 

(6.53 ft) 


2.1 Model Configuration and Instrumentation 

The one-seventh scale F-5A model was modified as necessary to correspond to 
the F-5A airplane on which dynamic pressure measurements have been made in flight. 
Specifically, removable wing-tip stores were mounted. Since local irregular- 
ities are detrimental to pressure measurements, the model surface was 
smoothed where necessary. Discontinuities greater than 0.02 mm (0.0008 inch) 
were eliminated, and surface waviness was minimized such that waviness height 
to waviness length did not exceed 0.002. A typical measurement of the wing 
top surface profile along Wing Station 8.00 (inches) is shown in Figure 2, 

A new leading edge flap at one deflection angle was provided on the 
right wing with instrumentation. This flap provided a smooth transition 
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into the wing contour and did not simulate the full-scale aircraft in 
having rearward and forward steps at this transition. (The same leading 
edge flap deflection was available for the left wing through use of an 
existing part. ) The new leading edge flap (5® deflection at root section) 
was attached to the wing by means of existing attach areas. 

The aileron can be deflected to various angles using existing brackets. 

The trailing edge flap is capable of being deflected in 5-degree increments 
down to 40® through use of a serrated disc. The left and right horizontal 
tails can be deflected from -20® to +10° in 5-degree increments through 
use of a serrated disc. Typical boundary layer transition strips were 
installed on the test model during the buffet test program. 

The F-5A scale model was equipped with 28 static pressure taps and an 
equal number of dynamic pressure transducers located closely to each other. 
Because of certain cutouts on the scale model, some transducers were located 
in slightly varied chordwise positions as compared to the airplane installa- 
tions. The static pressure taps consists of 0.76 mm inside diameter tubes 
terminating flush and square with the surfaces, and were connected to a 
six module scanivalve located at the front portion of the fuselage internal 
cavity. The dynamic transducers were Kulite Model LQ2-156W-4 with closed 
reference chambers. Specifically, 21 pressure taps and Kulite transducers were 
installed on the right wing, of which 18 taps and dynamic pressure transducers 
were located on the wing top surface with the ricmaining three sets located on 
the bottom surface. Three taps and dynamic transducers were located on the 
top surface of the left wing. Four taps and dynamic transducers were installed 
on the right horizontal tail with two sets each on the top and bottom surfaces. 
The transducers were mounted in such a way as to minimize model strain effects 
on the dynamic pressure measurements. For this purpose, the cutouts and wiring 
paths were designed to maintain the structural integrity of the model. Electri- 
cal leads from the transducers were connected to shielded-pair wires (with shields 
insulated from the model) as close to the transducers as practical while retain- 
ing good accuracy in the force measurements. The dynamic instrumentation was 
bonded flush with the wing surface with epoxy. Pressure taps were soft-soldered 
in the leading edge flap. No pressure taps were installed on the right wing 
leading edge flap with 0° deflection. 
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In addition, a pair of semi-conductor type bridge gages (Kulite 
ULP— 120— 160) were installed on both wing root sections to measure the bend- 
ing and torsion moments of the wing under dynamic loads. These were 
designated as lb, It, for the right wing and 2b, 2t, for the left wing. 

Three accelerometers (Endevco 2264—150) were installed in the model GG 
location and both wing-tips inside the sidewinder missiles. The accelero- 
meters were designated as la, 2a, 3a, with la located at the right wing-r- 
tip, 2a at the left wing-tip, and 3a at the CG, Transition strips were 
installed on the wing and tail surfaces of the scale model at approxi- 
mately 107 q chordwise positions. The detail instrumentation locations on 
the wing are illustrated in Figure 3, In the figure, a number inside 
the parentheses indicates that the corresponding transducers are on the lower 
lifting surface, while a number with no parentheses indicates that the trans- 
ducers are on the top surface. Table 1 gives the instrumentation information 
in a tabulated form. 

2.2 Model Support Systems Corresponding Vibration Modes 
Fixed Support System — Phase I 

As described previously, during Test Phase I, a standard sting was 

used to support the model. A six— component balance was installed inside 

the model fuselage in front of the sting mounting system. The outputs 

from the balance were designated as In, 2n, for the side forces, r for the 

roll moment and 0 for the axial force. After the F-5A scale model was 
a 

installed in the Eleven— Foot Transonic Tunnel, the model and the support 
system were subjected to sinusoidal excitations to determine the modal 
characteristics of the assembly. The modal data obtained during the test 
were recorded graphically in Figure 4« In addition to the natural modes 
of Figure 4, a support system bending mode was observed at 6.6 Hz. A 
support system torsion mode was observed at 105 Hz. In both cases, the 
model behaved as a rigid inertia mass. A mode corresponding to balance roll 
was observed at 13.6 Hz. 
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STATIC & DYNAMIC PRESSURE TRANS 
DUCERS (NUMBER IN PARENTHESES 
INDICATES TRANSDUCERS ON LOWER 
WING SURFACE) 
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Figure 3, Locations of Dynamic Instrumentation on F-5A Scale Model 



TABLE 1 F-5A SCALE MODEL INSTRUMENTATION 


RH WING INSTRUMENTATION 


SYMBOL 

WS+ 

1 

1 

Y+ 

STATIC 
PRESS. TAP 

PRESS. 

TRANSDUCER 

DIST. 

FROM LE 

IS 

1 

3.259 

(128.31) 

. .1102 

(4.34) 

2S 

2 

It 

fi 

.2646 

(10.416) 

3S 

3 

It 

If 

.4413 

(17.374) 

4S 

4 

H 

ff 

.6618 

(26.054) 

5S 

5 

If 

If 

.8273 

(32.571) 

6S 

6 

2.802 

(110.30) 

.3142 

(12.369) 

7S 

7 

If 

If 

.7140 

(28.112) 

8S 

8 

Tl 

If 

1.071 

(42.168) 

9S 

9 

2.362 

(93.00) 

.1741 

(6.853) 

lOS 

10 

If 

If 

.3656 

(14.392) 

IIS 

11 

11 

If 

1.044 

(41.111) 

12S 

12 

1.814 

(71.40) 

.9375 

(36.911) 

13S 

13 

If 

If 

1.278 

(50.323) 

14S 

14 

If 

fi 

1.704 

(67.102) 

15S 

15 

1,270 

(50.00) 

.2518 

(9.912) 

16S 

16 

II 

fi 

.5035 

(19,824) 

17S 

17 

II 

II 

1.057 

(41.622) 

18S 

18 

II 

II 

2.014 

(79.289) 

(19S) * 

(19) * 

3,259 

(128.31) 

.2646 

(10.416) 

(20S) * 

(20) * 

II 

If 

.6618 

(26.054) 

(21S) * 

(21) * 

2.362 

(93.00) 

1.044 

(41.111) 

BENDING 
STRAIN GAGE 

TORSION 
1 STRAIN GAGE 



lb 

i 

0.773 

(30.45) 

1.234 

(48.59) 

ACCELEROMETER 



la 


3.848 

(151.50) 

-.2934 

(-11.55) 


10 

24 

40 

60 

75 

22 

50 

75 

10 

21 

60 

44 

60 

80 

10 

20 

42 

82 

24 

60 

60 








tabu: 1 F-5A SCAI£. MODEL INSTRUMENTATION (CONTINUED) 


LH WING INSTRUMENTATION 


SYMBOL 

ws"*” 


Y+ 


STATIC 
PRESS. TAP 

PRESS. 

TRANSDUCER 

DIST, 

I 

FROM LE 

fo 

CHORD 

26S 

26 

3.259 

(128.31) 

. 2646 

(10,416) 

24 

27S 

27 

M 

IT 

.6618 

(26.054) 

60 

28S 

28 

2.362 

(93.00) 

1.044 

(41.111) 

60 

BENDING 
STRAIN GAGE 

TORSION 
STRAIN GAGE 




2b 



_ 

0.773 

(30.45) 

1.234 

(48.59) 

43 

ACCELEROMETER 




2a 

3.848 

(151.50) 

-.2934(-11.55) 

— 


RH HORIZONTAL TAIL INSTRUMENTATION 


SYMBOL 

WS^ 

Y+ 

DIST. FROM LE 

7o 

CHORD 

STATIC 
PRESS. TAP 

PRESS. 

TRANSDUCER 

22S 

22 

1.814 

(71.40) 

,2757 (10.85) 

40 

23S 

23 

1.270 

(50.00) 

.4368 (17.20) 

40 

(24S) * 

(24) * 

1.814 

(71.40) 

.2757 (10.85) 

40 

(25S) * 

(25) * 

1.270 

(50.00) 

.4368 (17.20) 

40 


* () Indicates the transducer mounted on the lower surface, 
t All dimensions are meters for the full scale aircraft. 

The corresponding values in the parentheses are inches. 
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Figure 4. Vibration Modes of F— 5A Scale Model Mounted on a Sting in Ames 
Eleven- Foot Transonic Wind Tunnel - Test Phase I 



Flexible Support System — Phase II 


The full-scale F-5A displays a wing- rock natural frequency of about 0,5 to 
1.0 Hz. This corresponds to a model frequency of 3.5 to 7,0 Hz. Since it is 
infeasible to allow the model to rotate freely about the roll axis in a 
wind tunnel test, a support system with a flexible roll constraint was 
designed for Phase II of the test program. This roll constraint was 
designed to have a natural frequency approximately equal to the natural wing- 
rock frequency so that impedance to the desired roll mode is minimized. The 
system allowed the model to oscillate +30° in roll. An adjustable damping device 
and a roll angle measuring potentiometer also formed part of the assembly. 

Except for the elimination of the balance unit, the pressure and response data 
instrumentations of the model remained unchanged from the first test phase. 
Dynamic data such as the roll angle, the model pitch and yaw oscillation angles, 
and the damping coefficient of the damper were recorded. 

The flexible support unit (Figure 5> was assembled with a specially designed 
hollow sting whose external dimensions were comparable to the standard 
sting used in the Eleven-Foot Transonic Tunnel. The overall length of 
the sting was 2.553 m (100.5 in.). At approximately one-third the length 
of the sting measured from the front end, a steel insert was fitted inside 
the sting which was used to anchor the rear end of a torsion bar. The 
torsion bar^designed to oscillate with the scale model along the roll 
axis with natural frequency of 3.5 Hz, had splines cut at both ends. 

With the rear end anchored to the sting insert as described above, the 
front end of the torsion bar extended approximately 0.1842m (7,25 in.) 
from the front end of the hollow sting. 

The front portion (0.324 m or 12.75 in.) of the sting was machined to 
accept two needle bearings and two lightweight thrust bearings which were 
fit into the inside of a sleeve shaped spacer. After the bearings were 
assembled, the spacer was fixed to the model fuselage body through a 1.27 cm 
(0.5 in. ) pin. 
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Beyond the bearing subassembly, the front tip of the sting was fitted 
with internal splines. The splines were used to accept the rear end of 
the damper subassembly shaft. The damper subassembly, measured approxi- 
mately 0.1651 m (6.5 in.), consisted essentially of the internal shaft, 
and the external housing and their attachments. The internal shaft was 
tied through the external splines to the sting front end. It was thus 
stationary along the roll axis. The external housing was fixed to the scale 
model through locating pins. The damper shaft was assembled with a 240® 
partial cylinder to form a 120® fan shaped cavity on top of the shaft. 
Correspondingly, the damper housing was attached with a fan shaped piston 
which covered an arc of 60®, When properly assembled with a number of 
seals and spacers, etc., and hydraulic liquid added to the cavity, the 
damper housing may rotate +30® about the shaft subassembly, while the 
hydraulic liquid was forced to move from one side of the piston to the 
other through three metering holes cut in the piston. The sizes of the 
metering holes which allowed the liquid to move were controlled by a 
spool shaft. By adjusting the lengthwise position of the spool shaft, 
the size of the metering holes were varied and thus the amount of damp- 
ing supplied by the damper unit could be controlled. The damper unit was 
locked out and the model prevented from roll oscillation when the metering 
holes were completely blocked by the spool shaft. The linear movement of 
the spool shaft was controlled by a Globe motor through a screw jack device. 
The fast acting lock-out mechanism of the spool shaft was activated by twp 
solenoids through a rocker arm. 

In front of the damper unit, and behind the two solenoids assembled 
inside the model fuselage, an angular potentiometer was used to measure 
the roll angle of the scale model relative to the sting. 

In addition to the angle-of-attack measuring device at the base of 
the sting, a danglometer was installed on top of the sting at a length- 
wise location about one-third from the front end of the sting. Strain 
gages were installed on the sting which were calibrated to determine the 
sting deflections in both the pitch and yaw planes. 



The calibration curves of the hydraulic damper with the 
scale model installed on the sting assembly obtained before and during 
the tests are shown in Figure 6 . The modal data of the assembly (the 
model and the support) during the test were recorded graphically in 
Figure 7, In addition, the first sting bending mode was observed at 
6.7 Hz and the nominal natural frequency of roll oscillation was 3.5 Hz. 

2,3 Test Points 

This program included the major test conditions conducted in the 
flight test program. These flight test conditions were M = 0,75, 0.85, 
and 0.925 for the wing leading and trailing edge flaps at 0® each and at 
4° and 12®, respectively. The wind tunnel program substituted a wing leading edge 
flap deflection of 5° for the deflection of 4° flight tested. The wind tunnel 
model did not provide a 4° leading edge flap deflection. It was expected 
that the difference of 1 ° in the deflection angle will not affect the test 
data significantly. 

The wind tunnel program also included a buildup of flap deflections. This 
was conducted to investigate the individual effects of each control surface 
and to provide the basis to measure the interaction effect between the two 
control surfaces. For Phase I testing, 52 test runs altogether (i.e., test 
points) were performed encompassing 11 model configurations defined by the con- 
trol surface deflections, etc. (leading edge flap 6 ^, trailing edge flap 5 ^, 
aileron 5 ^, horizontal tail surface 5 j^, wing-tip missiles on or off, wing 
transition strips on or off). Three Mach numbers were selected for testing 
(0.75, 0.85, 0.925) which were nominally identical to the numbers used 
previously in a flight test program. The nominal Reynolds numbers used in the 
first phase were 13.94 x 10^, 12.30 x 10^, 11.48 x 10^ and 7.38 x 10^ per m. 

Based on the model mean aerodynamic chord of 0.3366 m, the corresponding 
Reynolds numbers were 4,68 x 10^, 4.14 x 10^, 3.86 x 10^, and 2,48 x 10^. 

Another test run parameter was the sideslip angle, which was set to either 
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CURVE II 
CURVE III 
CURVE IV 


RUNS 1-16 
RUNS 17,18 
RUNS 19,20 
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(RUN 24 CALIBRATION DATA UNCERTAIN) 
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Figure 6, Calibration curves of the 
Hydraulic Damper with the 
F— 5A Scale Model Installed 
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Figure 7. Vibration Modes of F-5A Scale Model Mounted on a Flexible Support 
in Ames Eleven-Foot Transonic Wind Tunnel - Test Phase II 



••0° or +8° for each individual run. All tests were conducted as pitch polars, 
i.e., the angle-of-attack was varied while the sideslip is held constant. 
Corresponding to each test point, the angle-of-attack was stepped up from 0® 
through 4®, 6®, 8°, 10®, 12®, 14®, 16®, and the data were recorded for all <y's. 
A description of the Phase I test points is given in Figure 8. 

As described previously, the second tunnel entry (Phase II) was 
carried out with a flexible support system that allowed for model oscilla- 
tions, especially those degrees of freedom corresponding to the Dutch roll 
motion. A total of 29 runs were made, identified as test points 1-29. 

Except for the support flexibility, the Phase II run conditions were 
essentially duplicates of those of Phase I. A constant Reynolds number 
of 7.38 X 10^/m was used throughout Phase II. The Phase II test points 
and their corresponding test point numbers of Phase I are given in Figure 9. 
Furthermore, Table 2 describes the test points or pairs of test points used 
to investigate certain effects during parametric studies along with the 
corresponding subsection numbers as they appeared in this report. 
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Section 3 


PRESENTATION OF TEST DATA 

All the testing specified in Figures 8, 9, were performed at Ames 
Research Center Eleven-Foot Transonic Wind Tunnel. After a scan of test 
results, the recorded force and static pressure data were reduced to co- 
efficients and the dynamic pressure and response data were stored on the fre- 
quency-modulated magnetic tapes for subsequent real-time and spectral pro- 
cessing. In this section, the processing and evaluation of the real-time 
test data as well as the static and dynamic pressure distributions are 
presented, followed in a subsequent section by the spectral analysis of 
the test data. 

3.1 Static Pressure and Dynamic Pressure Distributions 

The static pressure and dynamic pressure data were examined first 
for correlation with buffet intensity and with onset of wing-rock. The 
tunnel data acquired during Run 8, Phase I, are presented below. Run 8 
data will be used as the base line data to evaluate the various effects 
to transonic buffet due to configuration and flight condition changes. 

Run 8 of Phase I test was conducted in test model configuration 2 

(Figure 8). This configuration included the tip-mounted Sidewinder 

missiles. The leading edge and trailing edge flap angles were 5®/12°. 

The horizontal tail surface and aileron settings were 0®. The nominal 

sideslip angle was 0°. The run was conducted at M = .926 with a Reynolds 

6 6 

number of 7.402 x 10 /m (2.256 x 10 /ft). The dynamic pressure was 
2 

18.232 KN/m (380.8 psf). The Reynolds number of the model, based on 

6 

the mean aerodynamic chord, was 2.49 x 10 . 

Run 8 covered the angle-of-attack settings from 0° to 16° in seven 
steps. Typical static and rms dynamic pressure coefficient distributions 
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on the right wing upper surface acquired during the run are presented in 

Figure 10 in this figure, the termination points C and C- at L.E. and 

T.E. were arbitrary and should not be taken literally. Referring to 
Figure 10, the development of the supersonic flow and the accompanying 
shock- induced separation can be seen in the variation of static pressure 

coefficient C along the chordwise direction with increasing angle"“of “attack, 
P ^ V 

The two conditions (supersonic flow and shock— induced separation^ corres- 
ponded to the high negative levels and subsequent flow recompression. 

The shock-induced separation moved upstream with increasing angle-of-attack. 
The increase in rms pressure coefficient C- in the shock/boundary layer 
interaction and downstream separated region can be clearly seen. 


It was noted that as the separation region developed over the upper 
wing surface with increasing angle-of-attack, the condition which 
corresponded to the separated flow developing over a "substantial*' part 
of the outer wing occurred over a relatively narrow angle-of-attack range 
(Aa ~ 2°). For tests conducted at a fixed angle-of-attack within this 
range, the dynamic pressures were found to be highly non-s tat ionary , 
particularly, in the shock boundary layer interaction region. This be- 
havior may be observed in the real time dynamic pressure traces as 
illustrated in the next subsection. 


3.2 Real Time Data for Dynamic Pressures and Model Responses 

The real time dynamic pressure traces of Run 8, Phase I, are 
presented in Figures (11-13) for transducers 2, 3, and 4 which are 
located at W.S. 128,31 corresponding to the 24, 40, and 607o chordwise 
stations, respectively. Figure 11 is for = 8° in which the shock 
interaction and the separated flow condition occurred in the wing trail- 
ing edge region as described previously (compare with Figure 10). The 
incipient shock interaction point was very close to transducer No, 4 as 
indicated by the random positive peaks in Figure 11. It was clear that this 
condition, a degree of non-s tationarity existed in the location of the 
shock interaction region. 


31 




Figure 10. Static Pressure (C^) and RMS Pressure C— ) Coefficient Distributions on the Right Wing 

Upper Surface, Run 8, Phase I - Configuration 2, p = O'”, M = .926, 6 /6^ = 5®/12® 

n f * 

6 = 0° (Continued) 
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Figure 11. Dynamic Pressures of Transducer Nos* 2, 3, 4 of the F5-A Model, 
Run 8, Phase I - Configuration 2, ~ 8°, P= 0°, M = .926, 


ir 
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Figure 12. Dynamic Pressures of Transducer Nos, 2, 3, 4 of the F-5A Scale Model 
Run 8, Phase I - Configuration 2, O' = 10°, p= 0°, M = .926, 
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By increasing the angle— of-at tack by 2® to O' = 10®, Figure 12 shows 
the shock interaction point had advanced to the 207o chord point (trans- 
ducer No. 2). However, the flow was highly random with the shock 
interaction region oscillating between the 20 and 407« chord points. 

Figure 13 shows the same data at an angle-of-attack (df = 16®) in which 
the outboard wing was fully stalled (compare with Figure 10). The 
turbulence can be seen to possess a more stationary character at this 
condition. In a larger time scale, the non- stationary fluctuation of the 
rms dynamic pressures for transducers 2, 3, 4, corresponding to = 10° 
is shown in Figure 14. The rms values of this figure were obtained by 
the integration of data within a traveling window 40 ms in width. As will 
be described later, the power spectral processing of the real time pressure 
data was applied to a finite and variable time span according to the 
frequency range of processing. In the high frequency range (up to 10 KHz), 
data from a relatively small time span was used for each ensemble. Aver- 
aging of the data from a number of ensembles yielded the power spectral 
density. A progressively longer time span was used for PSD processing in 
the lower frequency ranges. As a result, the non-stationari ty of the data 
may cause substantial deviation and uncertainty in the spectral density 
functions for conditions such as the case of = 10® depending on the frequency 
range and the specific ensemble time span selected for processing. 

In the second test phase, a special sting was designed incorpora- 
ting a torsional spring and damper which allowed the model to oscillate 
in roll at a natural frequency simulating the Dutch roll motion during 
wing-rock. The maximum roll angle was 30® single amplitude. In order to 
demonstrate the effect of rolling motion, typical response data of Run 5, 

Phase II, at a = 10° as recorded on an oscillograph are shown in Figures 15 
and 16. The test conditions are described below. 

Configuration 2 , 6 /6 = 5°/12° 

n f 

6 = 0 °/ 0 °, 6 ^ = 0 ° 

a h 

M = .925 

= 7.382 X 10^/m (2.25 x 10^/ft) 

2 

Pq = 29.80 KN/m (622.47 psf) 

Q = 18.20 KN/m^ (377. 70 psf ) 
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The nominal natural frequency of roll oscillation was 3.5 Hz, The damping 

ratio (c/c ) of the roll device was approximately 0.05. With increasing O', 
cr 

the model started to oscillate with substantial amplitude at a = 10°. 

Typical response data obtained at a = 10° as recorded on an oscillograph 
are shown in Figures 15 and 16. In Figure 15, the roll angle and the 
pitch angle oscillations due to sting bending are presented. Evidently 
the oscillations started in a random manner, reached a peak of approximately 
20° single amplitude, and then died out. The major frequency of oscillations 
was 4.0 Hz. The pitch angle oscillations had the same major frequency and 
were in phase with the roll oscillations. The random appearance of the 
roll oscillations and the subsequent reduction in oscillation amplitude 
were typical of the model motion at and above the critical angle-of-attack 
throughout the test phase. The response in this case was of the random type, 
as distinguished from the cyclic or periodic type (see, for example, Figure 39, 
page 77 of this report). For the case of random responses, the buffet pressures 
excited the model in a random fashion. When the narrow band pressure spec- 
tral density reached a certain level, the model was driven to oscillation of 
ever-increasing amplitude. After the oscillation reached a certain amplitude, 
either the buffet pressure intensity dropped, or the aerodynamic damping of 
the model reached such a level as to drastically reduce the oscillation ampli- 
tude as shown in Figure 15. The time of appearance and the length of the 
sustained oscillations were both random. On rare occasions, the model os- 
cillated continuously without letdown until the test condition changed. 

Also, the highest amplitude of oscillations occurred at or near the critical 
angle-of-attack The oscillation amplitude usually dropped with increasing 

angle-of-attack beyond cr’' . This behavior was observed even when the dynamic 
pressure PSD level continued to increase with increasing O'. Covering the 
identical time span of Figure 15, Figure 16 shows the right wing root 
bending moment, the right wing-tip rms acceleration and the CG acceleration 
time histories. An application of the Randomdec program to determine the system 
aerodynamic damping corresponding to high-amplitude roll motion will be presented 
later in this section. 

3.3 Comparison of Pressure Distributions Under Various Model 
Configurations and Test Conditions for the Fixed Model 

Evaluation of test data has shown that significant flow pattern develop- 
ment on the aircraft lifting surfaces and the resulting changes in the 
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Figure 15. Oscillographs of F-5A Model Roll Angle and 
Pitch Angle Recorded During Run 5, Phase II, 

Of *= 10\ p = o\ M = 0.925, 6. /6^ = 5°/12°, 

6^ = O'* (Timing marks appeared on bottom trace 
at one second Intervals.) 







static and dynamic pressure distributions occurred in the angle-of-a ttack 
range 6 degrees to 12 degrees. A detail study on the pressure distribution 
provided insight to the aircraft buffet mechanism. In the subject test 
program, 52 separate runs were made during the first test phase. Parametric 
changes in either the model configuration or test condition were carried out 
systematically. By comparing the pressure patterns corresponding to two 
test runs where the difference was limited to a single test parameter, the 
net effect of the changing parameter to aircraft buffet may be identified. 

In the following, comparisons were performed on a number of test run pairs 
based on the static and the dynamic pressure coefficient distributions. 
Effect of flap settings will be presented in the next section together 
with the power spectral density plots. 

Mach Number Effect 


In Figure 17 , the dynamic pressure coefficient distributions for 

Runs 19 and 22, Phase I, are presented. Runs 19 and 22 were conducted 

under model configuration 3, i.e., 6 /6 = 5°/12°, 6 = 0°, P = 0°, 6, = -10° 

n f a h 

Tip-mounted missiles were present in both runs. Run 19 was conducted at 
M = 0.925 and Run 22 at M = 0.75. The comparison of the two runs demon- 
strated the classical Mach number effect on aircraft wing flow separation. 

At M = 0.925 (Run 19), the shock boundary may be traced through the peak 
dynamic pressure locus. It was initiated near the trailing edge of the 
wing-tip region and gradually expanded toward the inboard and upstream 
region with increasing Of. The wing top surface was completely separated 
at or = 12°. For M = 0.75 (Run 22), the leading edge induced flow separa- 
tion appeared up to approximately 61 percent semi-span at = 6°. With 
increasing a, the flow of the wing top surface was completely separated. 
Throughout the process, no prominent dynamic pressure peaks were in evidence 
that identified the local shock. 

Sideslip Effect 

Figure 18 shows the dynamic pressure distributions for Run Nos. 8 
and 10. Both runs were conducted using model configuration 2, i.e., 

6^/S^ = 5°/12°, 6^ = 0°/0°, = 0°. The wing was attached with wing-tip- 
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mounted Sidewinder missiles. Both runs had a Mach number of 0.925 with a 
Reynolds number of 7.38 x 10^/4n (2.25 x 10^/ft). The sideslip angle ^ of 
Run 8 was zero, while in Run 10, p = 8°, the flow direction was somewhat 
closer to a line normal to the right wing leading edge. Furthermore, the 
vortex effect created by the front fuselage to the right wing was alleviated 
with a positive p. The two factors described above directed to a turbulence- 
free and more favorable pressure gradient along the chordwise direction, 
which tended to move the shock and shock— induced flow separation downstream. 
Throughout the a range presented in the figure, the shock and the shock- 
induced flow separation were retarded and pushed downstream corresponding 
to a positive p (Run 10). It was further noticed that beyond = 12°, the 
flow on the top wing was completely separated. As a result, the dynamic 
pressure distributions for the two cases became less and less distinguishable. 

Figure 19 is a set of similar pressure plots for Run Nos. 8 and 9. 

Run 9 was conducted with a sideslip angle P = “8°, which was expected to 
exhibit a reverse trend of Run 10 (p * 8°), i.e., forward movement of the 
shock front for the same a. Indeed, this was the case, especially for 
a = 10°. For O' = 6°, the almost complete flow separation along the wing 
chord at 85 percent semi-span was conspicuous for p = -8°. The exact 
reason for this behavior remained to be resolved. Additional discussion 
on the development of flow separation bubble on the scale model with side- 
slip will be presented later in this section. 


Effect of the Tip— Mounted Missiles 

Figure 20 presents the dynamic pressure distributions of Run Nos. 8 
and 27. Run 27 was identical to Run 8 except that the tip-mounted missiles 
were absent in the former case. An examination of the figure indicated 
that at a = 10°, the missiieless wing-tip caused the flow separation to be 
effectively suppressed. As soon as the suppressing effect was no longer 
dominating, the dynamic pressure patterns of the two runs then became almost 
indistinguishable (see Ck' = 12°). In other words, the tip-mounted missiles 
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had an adverse effect in promoting the flow separation at a certain critical 
a angle. This effect was greatly diminished after the britical angle-of- 
attack was surpassed. 

Tail Deflection Effect 

Both Run Nos. 39 and 49 were conducted with retracted leading edge 

and trailing edge flaps. The ailerons were also retracted. Both runs were 

made at M = 0.925, while Run 39 had a somewhat higher Reynolds number 
6 6 

(11.50 X 10 /m versus 7.38 x 10 /m for Run 49). The major difference between 
the two runs was that Run 39 featured a horizontal tail surface setting 
of 0 degree, while that of Run 49 was -10 degrees. 

In Figure 21, the dynamic pressures of the two runs are presented. 

In examining this figure, it was noted that the deflected horizontal tail 
tended to delay or alternate its shock-induced flow separation in that part 
of the wing directly upstream of the tail. Correspondingly, for a given 
angle-of— attack, the shock strength and the level of dynamic pressure behind 
the shock tended to be stronger in the outer half span of the wing. This 
effect was most prominent at a range 8-10®. Beyond this, the effect of the tail 
surface setting became less prominent. 

3.4 Comparison of Dynamic Pressure Coefficient Distributions 
For the Fixed and Flexibly-Supported Models 

In Phase I tests, the Northrop F-5A model was supported on a regular 
sting with a balance unit located inside the fuselage. In the Phase II 
tests, the same model was supported on a flexible system which allowed the 
model to oscillate along the roll axis under unsteady aerodynamic forces. 

A hydraulic damping device was installed to regulate the angle of oscilla- 
tion. The details of the flexible system are described in the previous 
section. 

In Figures 22 through 25, the dynamic pressure coefficient distri- 
butions on the right wing upper surface acquired during the two test phases 
under identical test conditions are presented. The plots cover the angle— 
of-attack range Of = 6®, 10°, 12°. The results of Phase I are plotted on 
the top of the figures with the corresponding data of Phase II plotted 
below. The test configuration and the run conditions are described in 
the following paragraph. 
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Figure 25. Comparison of Dynamic Pressure Coefficient Distributions for the 
Fixed Model (top) and the Oscillating Model - Configuration 2 , 





Figures 22, 23, 24: Configuration 2, 6^/6f = 5°/12®, 6h ~ “ 0®/0®, 

Sidewinder missiles at wing- tips, M = 0.925, = 2,485 x. 10®. Figure 25: 

Same conditions as above except that M = .75, 

Additional dynamic pressure coefficient comparisons for the fixed and 
flexibly-supported model were made in Figures 26 to 30 for the following cases: 

Figure 26: Configuration 3, 6^/ 6f = 5®/12“, 6]^ = -10°, 6a = 0°/0°, 

Sidewinder missiles at wing-tips, M = 0.925, = 2.485 x 10®, 

Figure 27: Configuration 8, 6^/6^ = 0°/0°, ~ 0°/0°, 

no Sidewinder missiles, M = 0,925, R^ = 2,485 x 10®, 

Figures 28-30: Configuration 9, same as Configuration 8 except that 

Sidewinder missiles were installed at wing-tips, M = 0,925, 
R^ = 2.485 X lO! 

The specific test phase and run numbers presented in Figures 22-30 are 
tabulated below: 



Test 

Phase-Run No. 

Test Phase -Run No. 

Side sli p Angle 

a* 

Figure 

22 

1-8 

II-5 

(3 = 0° 

10" 

Figure 

23 

I- 10 

1 1- 6 

P = 8** 

12" 

Figure 

24 

1-9 

II-8 

P = -8° 

11° 

Figure 

25 

1-17 

11-20 

P - 8" 

14" 

Figure 

26 

1-19 

11-14 

3 = 0° 

10° 

Figure 

27 

1-36 

11-29 

3 = 0° 

8° - 10° 
intermitten 

Figure 

28 

1-39 

11-21 

p = 0° 

8° 

Figure 

29 

1-41 

11-22 

3 = 8° 

10° 

Figure 

30 

1-40 

11-23 

P = -8° 

10° 


= the angle-of-ra ttack at which the flexibly-supported model started 
to oscillate along the roll axis. 

In examining the pressure distributions for the fixed model and the 
flexibly-supported model, the following conclusions may be drawn: 
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1. For angles-of-attack under the critical angle q^', it was expected 
that the dynamic pressure coefficient distributions would be identical in 
phases I and II under otherwise identical test conditions. The plotted data 
showed that corresponding to the two phases, the trends of the dynamic pressure 
coefficient distributions were the same which reflected the shock development 
(for M = 0^925) and the expansion of the flow separation region, etc. 

However, substantial deviation in local pressure did exist. For instance, 

in Figure 24, at o' == 6 degrees, the rms pressure levels along 61 percent and 
72,8 percent semi-spans were substantially higher for Phase II, Contrarily, 
in Figure 25, o = 10 degrees, the very high rms pressure peak at 20 percent 
chord and 72,8 percent semi-span location shown in Phase 1 was not observed 
in Phase II, These variations of the dynamic pressure data indicated a 
certain degree of randomness of the buffet phenomenon, with a predictable 
trend development for any given test condition. 

2. In general, at the critical angle-of-a ttack or' when the model started to 

oscillate along the roll axis, the model motion tended to make the pressure dis- 
tribution more even, spacewise. This condition was observed in a number of cases, 
but with prominent exceptions. One exception was Figure 23, = 12°, where the 

high dynamic pressure for the model with roll freedom may be contributed to the 
coincident location of the shock region and the pressure transducer( s) . 

3. For a > or' y the differences in the dynamic pressures of the fixed and 
oscillating models tended to diminish as the angle-of-a ttack was further 
increased. 

Referring to Figure 26, the extremely high C- value (0.2788) of transducer 
No. 1 at O' = 10 degrees indicated that the data was recorded prior to the large 
amplitude model oscillations, which also started at oT'' = 10 degrees. It can be 
visualized that an intense shock was originated at a location near transducer 
No. 1, which served as an anchoring point of the shock line on the wing 
surface. Toward the inboard stations, the lesser rms peaks indicated either 
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Figure 30. Comparison of Dynamic Pressure Coefficient Distributions for the Fixed Model (top) and the 
Oscillating Model - Configuration 9, P = -8®, M = .925, G /o^ = 0°/0°, 6 = 0°, 



reduced shock strength, or unstable shock location, or most possibly that the 
shock was not located precisely over the pressure transducer. Similar phen- 
omenon may be observed in Phase I, Run 39 data at a = 8 degrees (Figure 28) 
and Phase I, Run 41 data at a = 10 degrees (Figure 29). 

Referring to Figure 27 , corresponding to the model tested without tip- 
mounted missiles, the substantial differences in the shock line location between 
Phases I and II at 0( = 6 degrees indicated the sensitivity of shock location 
to minor variations in test conditions. The lack of the tip— mounted missiles, 
and their trailing vortex flow, seemed to aggravate the situation and make 
the shock location and the resulting flow pattern less stable, as compared 
to the cases with tip-mounted missiles. 

Comparison of power spectral density data for fixed and flexibly- supported 
models (Run 8 , Phase I vs Run 5, Phase II) will be presented in next section. 

3.5 The Development of Flow Separation Bubble 
on the Scale Model with Sideslip 

For the purpose of detecting a special trend in separation bubble develop- 
ment on the wing surface, a number of rms pressure coefficient distribution 

(C-) data were examined. Attention was directed toward the pressure data on 
P 

the right wing upper surface when the model had a positive sideslip (p - 8 
degrees). For those runs, the right wing was on the windward side. The runs 
under study and the test conditions are tabulated below: 


Run No. 

Run 10 
1 Phase I 

Run 28 
Phase I 

Run 38 
Phase I 

Run 6 
Phase II 

Configuration 

2 

6 

8 

2 

«£ 

50/120 

5®/12° 

0 ®/ 0 ® 

5®/ 0° 

«a 

0 ®/ 0 ° 

0 ®/ 0 ® 

0 ®/ 0 ® 

0 ®/ 0 ® 

«h 

0 “ 

0 ® 

0 ® 

0 ® 

6 

8 ® 

8 ® 

8 ® 

8 ® 

M 

.929 

.927 

.927 

.926 

Rj^, /m 

7.441 X 10^ 

7.379 X 10^ 

11.51 X 10 

7.382 X 10 

P^, KN/m^ 

30.35 

30.21 

47.83 

29.90 

P^, KN/m^ 

52.96 

52.62 

83.37 

52.05 

Q, KN/m2 j 

18.32 

18.17 

28.80 

17.96 
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In the four runs listed above. Run 10, Phase I and Run 6, Phase II were 

conducted with the wing-tip missiles installed. The tip missiles were off 

for the remaining two runs. The rms pressure coefficient distributions (C— ) 

for the four runs with variable a are plotted in Figures 31 to 34. It may be 

pointed out that C— is identical to <p>/q as defined by Mabey in Reference 4. 

In dealing with the general cases of the separated flow, Mabey noted that 

close to the reattachment point of a separation bubble, a maximum value of 

<p>/q or C— between 0.04 and 0.10 was usually recorded. For the F-5A data 

being processed, C— was in the order of 0.01 or less in the unseparated flow 

region. When the local flow became separated, a substantial increase of C— 

P 

value was always observed. The C— value in the separated flow region of the 

wing may vary within the range of 0.04 to 0.16. In isolated cases, the 

local C— value was as high as 0.28 which identified a strong local shock, 

see, for instance. Figure 26, a = 10°. On the horizontal tail, the maximum 

C— value was substantially higher than 0.16. It was found convenient to use 

the above-mentioned C— values recorded on the wing surface to trace the 

P 

development of a separation bubble and the region of the separated flow. 

The validity of this approach was confirmed by comparison of the flow sep- 
aration data developed in a F-5A flight test program using tufts on the wing 
surface . 

In Figure 31, five plots of the C— coefficients for Run 10, Phase II 

are presented corresponding to a = 4 degrees, 6 degrees, 8 degrees, 10 degrees, 

12 degrees. At C3t = 4 degrees, the single high value of C— at 72.8 percent 

P 

semi-span and 75 percent chordwise position (transducer No. 8, see Figure 3) 
indicated the appearance of a shock-induced separation bubble at that loca- 
tion. Judging from the proximity of the location to the trailing edge, and 
considering the subsequent flow developments, the separation bubble seemed 
to extend toward the trailing edge and no reattachment took place on the wing 
surface. The separation bubble remained at the same general location at 
a = 6 degrees but the C— value was further increased to 0.079. At a = 8 degrees, 
the separated flow region was extended inboard toward the aileron and the 
trailing edge flap, with additional chordwise expansion toward the leading 
edge at a = 10 degrees. Up to a = 10 degrees, the flow at 85 percent semi-span 
seemed to remain completely unseparated. When a reached 12 degrees, a 
major expansion of the separation bubble took place. Except for a limited 
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Figure 32. Dynamic Pressure Coefficient Distributions of the F-5A Scale Model, Run 28, 

Phase I - Configuration 6, 9 = 8®, M = .927, 6 /6. = 5Vl2% 6 = Missiles Off 







area on the extended leading edge flap, the flow on the top surface o£ the 
right wing became completely separated. In visualizing the development of 
the separated flow on the windward side of the semi-span in a sideslip con- 
dition, it may be noted that the slanted air stream flow direction (positive $ 
for the right wing) seemed to create a favorable flow condition in its wing-tip 
region to retard the local flow separation. Similar behavior may be noted, with 
individual variations, in the other three runs where the sideslip angle was also 
8 degrees (Figures 32 to 34). For Phase I, Run 28 (Figure 32), the separated flow 
development was very similar to Phase I, Run 10 described above. For Phase I, 

Run 38 (Figure 33), both leading edge and trailing edge flaps were retracted, the 
location of the initial flow separation seemed to be closer to the wing-tip 
as reflected in the C- distribution at gjp =4 degrees. The location of the 
maximum C- coefficient moved spanwise alternatively inboard and outboard as 
the angle-of-attack was increased. In comparing with Figure 32, where the 
leading and trailing edge flaps were extended (6^/fi ^=5° /12° ) , it seemed that 
the extended trailing edge flaps tended to keep the shock induced flow 
separation region toward the inboard portion of the semi-span at a moderate 
angle-of-attack (or = 8°). This can be attributed to the increased camber in the 
inboard portion for the extended flap cases. The initial development of flow 
separation for Phase II, Run 6 (Figure 34) was very similar to those of Phase I, 
Run 10 (Figure 31). This was as it should be since both runs featured the same 
model configuration except for the support flexibility. 


The dynamic pressure data for the leeward side of the wing semispan for 
the model with a sideslip were examined. In all cases, flow separation 
seemed to originate at the wing-tip region. This was in contrast to the wind- 
ward side, where flow separation was developed initially at a location other 
than the wing— tip region. Correlation study of pressure data under sideslip 
conditions will be presented in the next section. 


3.6 Determination of the System Aerodynamic Damping 
Using the Randomdec Program 

In the tests carried out during Phase II, as described previously in 
subsection 3.2, the major model motion was due to support flexibility which 
allowed substantial Dutch Roll type oscillations known to be important in 
aircraft buffet. As the angle-of-attack setting approached the critical 
angle-of-attack a*, the model started to oscillate with substantial amplitudes 
along the roll axis, with correspondingly lesser amplitudes in yaw and side- 
slip. The major frequency of the oscillation was approximately 4.0 Hz. On 
rare occasions, the model oscillated steadily and continuously through- 
out the recording period, such as Run 26 at ot = a* = 10°. The real time 
data of pitch, yaw, and roll oscillations are given in Figure 35. Under most 
test conditions, the motion was random. The oscillations started 
in a random manner, reached a maximum amplitude and then died out. Typical 
real time data of roll and yaw oscillations, obtained in Run 5, Phase II, 
at a = a* = 10° are shown in Figure 36a. In Figure 36b, the data in a shorter 
time segment (between the two arrowheads shown in Figure 35a) are repeated 
with an enlarged time scale. Based on the real time data, it was clear that 
the yaw angle oscillations had the same major frequency with a slight phase 
shift from the roll oscillations. The yaw oscillation amplitude was com- 
paratively low. It was believed that the aerodynamic damping of the model 
changed abruptly at the critical angle-of-attack a* which affected the os- 
cillation amplitudes. The Randomdec program was chosen to determine the system 
aerodynamic damping based on the measured roll oscillations. 

The Randomdec program has been used extensively in determining the 
damping of a major response mode under random excitation such as that ex- 
pected in transonic wind tunnel testing. For multiple-mode responses, proper 
filter techniques are used to enhance the processing of the primary mode 
data. To obtain a Randomdec signature, a number of data segments are collected 
each having the same initial amplitude (or other characteristics, such as 
the normalized initial peak, etc.). Ensemble averaging is performed on the 
accumulated data. If the system is linear and the excitation random and 
stationary, the ensemble average then converges to the transient response 
of the system, due to the selected set of initial conditions. 
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Figure 36a. Time Histories of F-5A Model Yaw Angle and Roll Angle Recorded During Run 5, 
Phase II, a = 10% 6 = 0°, M = 0.925, 6/6^ = 5°/12% 5 = 0°. (Timing marks 

appeared on top trace at one second intervals) 



Two sets of measured roll data were digitized to provide the input to 
the Randomdec program. For Run 26, a = a*= 10°, the model oscillated 
steadily and continuously throughout the recording period. For Run 5, 
a = a* = 10°, the motion was highly random. The measured roll data 
were digitized covering the complete recording period with a sampling rate 
of 521 samples/sec. A low pass analog filter with a cutoff frequency of 
220 Hz was used during the digitizing process. 

Run 26 of Phase II test was conducted in test model configuration 11, 

'5 /6 = 0°/0°, 6 = -10°, 6 = 0°/0°, 3= 8°, Sidewinder missiles at wing-tips, 

nr n 3 . 

M = 0.925, and a* = 10°. Since the steady roll motion of Run 26 had only 
one major frequency of oscillations at approximately 4.0 Hz, no special 
filter technique was required to obtain the Randomdec signature. Figure 
37 shows the Randomdec signature obtained using a data span of 34 seconds 
(or 17,775 samples) with 140 ensemble averages. The damping ratio (structural 
damping plus aerodynamic damping) was 0.000816. The natural frequency was 
3.976 Hz. Based on the calibration curves of the hydraulic damper with the 
F-5A scale model installed (Figure 6), the damping ratio of the model with 
air off is 0,012, As a result, the aerodynamic damping ratio for Run 26 at 
O' = 10® was estimated to be -,0112, 

For Run 5 at a = cy* = 10®, the lengths of the selected time segments for 
various phases (the initial, high-amplitude, termination phases) were rather 
short and the statistical errors are expected to be very high. For most time 
segments, difficulties were encountered in obtaining a meaningful Randomdec 
signature which might represent a free vibration decay curve even though the 
roll signal was bandpass filtered before determining the Randomdec signature. 

It was believed to be caused by the short time span used. However, a meaning- 
ful signature was acquired for one specific segment (as indicated by letter A in 
Figure 36a) with low roll amplitudes. The result is given in Figure 38, For 
segment A, a total number of 1626 samples (3,1 seconds) were used and 11 
ensemble averages were made. The natural frequency obtained was 3,48 Hz 
and the damping ratio was 0,0176, The corresponding structural damping ratio 
was 0,0575 for Run 5 and the aerodynamic damping ratio was thus -,0399. A 
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Figure 37. Randomdec Signature of Roll Oscillations of the F-5A Scale Model, 

Run 26, Phase II, a = 10®, 3 - 8“, M « 0.925, 6 /6. « OVO®, 6, = -10 



Figure 38. Randomdec Signature of Roll Oscillations of the F-5A Scale Model, 
Run 5, Phase II, a = 10°, 6 = 0°, M = 0.925, 6 /6 = 5°/12% 5, = 



bandpass digital filter was applied to the roll signal before determining 
the Randomdec signature. The cutoff frequencies of the bandpass filter were 
2.5 Hz and 3.8 Hz and the difference between a cutoff frequency and the 
frequency where the filter has attenuated approximately 6 dB was 0.8 Hz. 

The above work experience indicated that the Randomdec technique for damping 
evaluation is applicable to test cases with a stable primary mode of excitation. 
For the case of purely random excitation of the scale model under buffet, 
the success of the technique is dependent on the relative amplitude of the 
primary mode of excitation, the frequency of occurrences of the major exci- 
tation mode, the data time span, and maybe other contributing factors. 


3.7 Aerodynamic Hysteresis and the Limit Cycle Concept 

In Phase II, two types of model oscillations along the roll axis were 
observed. In one type, the oscillation was approximately a constant-amplitude 
sinusoidal motion. The corresponding local pressure data also featured 
a periodic pattern, even though highly nonlinear with respect to the model 
motion. In a second type, the oscillation was essentially sinusoi- 
dal, but the amplitude of the oscillation changed from cycle to cycle. The 
roll motion existed in a pattern of irregular bursts, and Interspersed by 
periods of near zero oscillation. The first type of oscillation was named 
the cyclic type; a typical case was Run 26, Phase II (illustrated in Fig- 
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(Figure 3i>) , with the following test conditions; 

6 /6^ = 0V0% 6, = -10° 
n I n 

6 = 0 °/ 0 ° 
a 

M = 0.925, Q = 17.95 KN/m^ (375.0 psf) 

= 7.382 X 10^/m (2.25 x 10^/ft) , or 2.48 x 10^ based 
on the mean aerodynamic chord of the model 

a = a* = 10°, 3 = 8°, 

The corresponding run in Phase I is Run 51, which had identical model 
configuration as Run 26, Phasa II shown above. It is noted that the Phase 
I runs were made with a standard sting with minimum roll flexibility, as 
against the case of Phase II where a special sting with substantial roll 
flexibility was used. The Reynolds number of Run 51, Phase I was 7.38 

X 10^. 

The second type of oscillation was named the random response type. A 
typical case was Run 5, Phase II^ (illustrated in Figure 36), with the test 
conditions as shown below: 

V^^f = ^"/12°, 6 ^ = 0° 

6 = 0 °/ 0 ° 
a 

M = 0.925, Q = 18.20 KN/m^ (377. 7 psf) 

= 7.382 X 10^/m (2.25 x 10^/ft), or 2.48 x 10^ 
based on the mean aerodynamic chord of the model 

a = a* = 10°, 3 = 0°, c/c^j.= 0.057 

The corresponding Phase I is Run 8, of which the Reynolds number was 7.38 x 10^. 

The real time roll oscillation and the corresponding fluctuating 
pressure data are further illustrated in Figure 39 for Run 26, Phase II. In Figure 
39, P^, P^ , P^ represent the local fluctuating pressures at 85% semi-span 
of the right wing top surface, with transducers located at 40%, 60% and 
75% chordwise positions. On top of the oscillograph plots, the local 
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Figure 39, Real Time Data of the Roll Oscillations and the Corresponding Local 
Pressures of Transducer Numbers 3, 4, 5 of the F-5A Scale Model, 

Run 26, Phase II, a = 10°, 6 = 8“, M = 0.925, 5^/6^ * 0°/0°, \ = - 
Q =» 17.95 KN/m^, c/c = 0.012. 






Aa oscillation was presented based on the following: 

Aa = - cos ot^ sin 3 + (^y/U) cos 3 (1) 

where <p is the roll oscillation angle, y is the spanwise distance from the roll 
axis, and U is the tunnel air speed. In examining the fluctuating pressure 
data vs, the local Aa changes, it was apparent that the sudden pressure 
changes, representing the loss of lift or lift recovery, and the high 
amplitude oscillations, representing the unstable shock movement during 
transition from local unseparated flow to flow separation (or vice versa) , 
followed a cyclic pattern which was synchronized with the local 
Aa changes. Specifically, when the local Aa was changed from negative to 

positive, the shock oscillations and the subsequent flow separation tended 
to lead the Aa change. When the Aa was changed from positive to negative, 
the local pressure change (a lift recovery) was usually more clear cut 
reflecting the momentum of the local flow moving the shock front rearward 
toward the trailing edge. The detail behavior of the fluctuating pressures 
and their impact on aircraft wing rock can be best studied together with the 
static pressure development on the wing surface. 

In Figures 40, 41, the chordwise static pressure distribution for Run 
51, Phase I and Run 26, Phase II are plotted side by side to illustrate 
the effect of roll oscillations to local static pressure distributions. 

Both figures are for 85% semi-span of the right wing. Figure 40 shows the 
upper surface pressure data corresponding to the nominal angle-of-attack 
range from 0® to 12°. Figure 41 shows the lower static surface pressure 
data for selected a’s when substantial changes in upper surface pressures 
took place. 

Referring to Figure 40, for Phase I corresponding to a standard rigid 
sting, the high pressure gradient near the midchord for a = 6°, 8° indicated 
a shock induced flow separation in the rear semi-chord region. At a = 10°, 
the shock moved forward and became less distinct. As a result, the pressure 
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Figure 41. Chordwise Static Pressure Distribution at the Bottom Surface of 85 Percent 
Semi span, F-5A Right Wing. The Left Plot is Based on Run 51, Phase I, Data 
The Right Plot is Based on Run 26, Phase 11, Mean Value Data. 



gradient was more evenly distributed along the chordwise location. With 
increasing angle-of-attack, a = 12® reflected a major breakdown of the un- 
separated flow. At this a, the leading edge induced flow separation caused 
a substantial loss in lift in the outer semi-span of the wing. (The actual 
flow breakdown happened at an a slightly greater than 10°.) 

The Phase II data of Figure 40 show the corresponding static pressure 
distribution at 85% semi-span when the model was free to oscillate along the 
roll axis. The test results indicated the model was near stationary up to 
ot = 9.6°. As to be expected, the static pressure data up to this a 
closely resembled the Phase I data. At a = a* - 10.3°, sustained roll 
oscillation took place with a roll frequency of approximately 4 Hz. 

Selected real time response and fluctuating pressure data at this a were 
shown previously in Figure 39. The ranges GL, G'L' on Figure 40 indicate 
the extent of the fluctuating pressure changes based on Figure 39 data. 

Beyond a - 10.3°, the roll oscillations had either reduced amplitude, or in 
some cases, completely disappeared. Correspondingly, the static pressure distri- 
bution at 85% semi— span followed a similar pattern to that acquired in Phase I. 

Corresponding to the critical a range when drastic pressure re- 
distribution took place on the upper surface of the right wing, the lower 
surface static pressure data are plotted in Figure 41 for both phases. In 
examining these data, and the real time fluctuating pressure data of the 
lower surface, it can be concluded that the lower surface flow was not 
affected to any appreciable degree by the model oscillation. In other 
words, no or little interaction between the lower surface pressure distri- 
bution and the aircraft oscillatory movements was observed. 

In order to investigate the static and dynamic pressure behavior on 
the left wing (the leeward side for 3 == 8°) corresponding to the same test 
condition (Run 26, Phase II) when sustained roll oscillations took place, 
the data of Run 27, Phase II for the right wing were examined (Figure 42) . 

Run 27 was performed with identical model and tunnel conditions of Run 26 
with the exception that 3 = -8°. The use of the leeward side (right) wing 
pressure data of 3 = -8° to replace the leeward side (left) wing pressure 
data of 3 = 8° was necessary because of the scant number of pressure trans- 
ducers installed on the left wing. The above substitution implied the 
assumption of perfect symmetry of the test model and the wind tunnel along 
their central sections. 
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Figure 42. Chordwise Static Pressure Distribution at Top Surface 
of 85 Percent Semispan, F-5A Right Wing Based on 
Run 27, Phase II, Mean Value Data, Q = 17,95 KN/m 


Referring to Figure 42, and comparing the data with the windward 
side data of Figure 40, it was seen that prior to complete- flow separation, 
the lift distribution on the leeward side was lower than that on the wind- 
ward side. This behavior partially contributed to the positive side- 
slip angle of Run 26 which acted with the wing sweepback to reduce the mean 
effective angle-of-attack on the leeward (left) side. At higher a, the 
vortex shedding and wake turbulence effect of the fuselage, and possibly 
other factors, contributed to the asymmetrical pressure distribution of the 
wing. In general, complete flow separation appeared earlier on the leeward 
wing in the upper surface near the wing tip. It was also noted that the static 
pressure on the leeward wing surface varied in a moderate manner during 
the roll oscillations at the critical angle-of-attack, a* = 10.3®. The 
above information indicated that the major contributing factor to wing rock 
of F-5A with sideslip was due to the pressure fluctuation on the top surface of 
the windward semi-span. 

Combining the mean static pressure data of Figure 40 corresponding to 
a = 10.3® for Run 26, Phase II, and the extreme fluctuating pressure data 
taken from points such as G, L as identified in Figure 39, a plot of the 
range of the chordwise pressure distribution during a roll oscillation cycle 
is shown in Figure 43. 

An examination of Figure 43 indicated dynamic overshoot of the static 
pressures for both positive and negative Acx. For negative Aa, the maximum 
lift was generated along the chord section higher than those recorded 
in steady state tests. The maximum lift occurred prior to the time when 
the minimum Aa value was reached, indicating a phase lead of 30® to 45®. 

Using the convention of Reference 8 for a overshoot, the following may be 
stated for the wing strip in the neighborhood of 85% semi-span. 

Aa = K [-<p cos a sin 3 + (^Py/U) cos 3] (2) 

s a o 

K = 2 
a 

where is the dynamic overshoot coefficient. The above equation and the 
pressure data as illustrated in Figures 39 to 43 indicated that the oscillating 
dynamic loads, if any, were dependent on the amplitude and phase of the 
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roll oscillations during wing rock. The dynamic loads and the motion thus 
formed an interacting pair. The existence of limit cycle type oscillations 
was determined by the energy balance of (1) the input to the aircraft in 
the form of fluctuating loads caused by mot ion- induced flow separation, 

(2) the energy consumed by aerodynamic damping when the aircraft oscillated 
in roll, possibly with other coupled motions. The energy input to the 
aircraft, essentially due to the fluctuating pressure loads on the top 
surface of the right wing for 3=8°, may be computed based on the Cp 
data vs the local translatory motion due to the roll oscillation. In 
Figure 44, the local data for a typical station are plotted vs the 
translatory motion h and the local angle of attack ot + Aa. To estimate 
the work done by the fluctuating forces in one oscillation cycle, the 
following assumptions are made: 

(1) The work was essentially in the form of the fluctuating pressure applied 
to the top surface of the windward semi-span, and may be represented by 

the loads on a wing strip with width Ab and located at 85% semi-span. 

, , , . ^ ^ , , , i(jJt iajt 

The corresponding local oscillatory motion is he = <pye 

(2) The Cp contour along the chordwise direction was defined by the upper 
and lower limit values of Figure 43. The fluctuating pressures varied 
in coordination with each other so that the limit cycle loop as 
described in Figure 44 for was followed by all pressures along 
the chordwise direction. Admittedly, this is a simplifiec assumption 
which served only to describe the fluctuating pressure trend. 

Applying the above assumptions, the work input per oscillation cycle 
simulating wing rock is; 


Wk 



Ab A Q dh 


( 3 ) 


A 




( 4 ) 
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where Q is the dynamic pressure, 0,0^ are the lower and upoer values 
as defined in Figure 43. A is the area covered by and and integrated 
from the leading edge to the trailing edge, using the^scale oS the right 
hand of the plot. The Integral ^^covers the complete oscillation cycle. 

Based on the scale relations between the model and the aircraft, if the angle 
of oscillation is the same for the aircraft and the model, the work input 
per oscillation cycle as shown in equation (3) may be converted to aircraft 
scale as follows: 


Wk . ^ = L X (Q . j ,) X Wk , - 

aircraft aircraft model model 

where L = 7 is the model linear scale factor. In the following, all the data 

are referred to the scale model results. 


The energy consumed by aerodynamic damping, and in the case of the test 
model, the additional mechanical damping built into the support system, are: 

(5) 

S ■ % <5 S b^/2U 

where is the modal mass in roll, g^ is the mechanical damping coefficient, 
is the circular frequency, <P is the roll oscillation amplitude, C^p is 
the roll damping coefficient, S is the wing surface area, b is the aircraft 
wing span and U is the free-stream velocity. The energy balance condition 
for constant amplitude wing rock oscillation is then: 



Ab AQ dh =^Xgj^ 


- L ) ^ d<p 
P 


( 6 ) 


In the above equation, the value Ab, or the effective width of the wing tip 
area with significant fluctuating pressures to supply work to the limit 
cycle, is treated as a parameter. This parameter is determined by the energy 
balance condition. The numerical values of the parameters used to calculate 
the energy balance condition for Run 26, Phase II, a * a* “ 10® are shown below. 


- 0.00437 m^-kg , “ 2 c/c^^ = 0.024 
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=■ 25.13 rad/ sec 


Q = 17.95 KN/m^ 


S » 0.3244 , b = 1.1 m 

U » 299 m/sec , c- - -0.135 

Xrp 

dh =* 0.0163 m 

^2 -1 

’'’’nax “r = 2.17 sec 

Based on the above data, the following estimate of energy balance is reached: 


0.1076 m 


<p =9 
max 


I' 


E * 6.159 tn-N/cycle 

(of which 43.5% Is due to mechanical damping, 56.5% is due to aerodynamic dampln 
Wk « 0.3445 Ab ra-N/cycle 
Ab = 0.35 X (|) - 0.1923 m 


The real time roll oscillation and the fluctuating pressure data of 
Run 5, Phase II are illustrated in Figure 45. The approximate Aa data 
due to the local downwash and upwash strokes are generated and plotted 
on the top of the figure. P^, P^, P^, P^ were located at 85% semispan 
of the right wing top surface, and at 10%, 40%, 60%, and 75% chordwise 
positions respectively. Run 5 was conducted with zero 6. The roll 
oscillation amplitude was substantially lower than the corresponding value 
of Run 26, Phase II (see Figure 39), partially because of the higher mech- 
anical damping set for the support system. The oscillations were irregular 
and intermittent. In examining Figure 45 data, it was noted that the flow 
at was unseparated for the majority of the time, with intermittent lift 
losses which were synchronized with ^ changes. The continuous pressure 
fluctuations at P^, P^, P^ indicated a separated flow mostly induced by 
a shock located somewhere between P^ and P^. The lift losses at P^^, i.e., 
sudden pressure increases, appeared when Aa = 0 and ot = +, a condition 
inducing leading edge stall and local flow separation. The lift losses at 
Aa s 0 indicated a phase lead of the flow separation caused by the dynamic 
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effect. The lift recovery at appeared at even though Aa may still 

be positive at the time. The overall picture was a state of moderate in- 
stability, with the shock line located close to the leading edge, 4n inter- 
mittent movement of the shock l^-ue toward the leading edge caused the complete 
stall of the local section. The local stall was partially recovered within 
one half oscillation cycle or less. Typical limit cycle data for are 
presented in Figure 46. In the figure, the local Cp vs. data are also 
plotted. The very limited changes restricted the C variation away from 
the ramp representing the steady state loss of lift. The narrow Aa range 
for the test case was considered restrictive, not being a typical free 
flight condition. In a free flight, any roll oscillations would introduce 
coupled yaw and sideslip oscillations, which in turn would cause substantial 
Increase in Aa range when the nominal angle-of-attack ot^ was high. 

It was also noted that if the random turbulences were ignored, the cyclic pres- 
sure fluctuations (being the source of the limit cycle ener<^y input) , appeared 


only near the outer span leading edge area. This was in co..crast to Run 26, 

Phase II case where the pressure fluctuations appeared in substantial amnlj.- 
tude throughout the majority of the chord length (see Figure 39) . As a 

result, the wing rock type oscillations of Run 5, Phase II, representative of cases 
of zero e for the test phase, were intermittent and limited in oscillation amnlitudc 


The energy balanced condition for a cycle of nonstationary wing rock 
oscillation such as cycle A and B of Figure 45 may be written as: 


Mk - i Mr “r - Lp) V d<P (7) 

where TJk is the work input per oscillation cycle, <p^ are the initial 

roll oscillation amplitudes (at ^ = 0) of the oscillation cycle. The first term 
Iti the right hand side of the above equation denotes the potential energy 
change of the flexible roll device. Similar assumptions may be made for the 
work input as in the stationary (3 ^ 0) case (see Equation 3), with the 
exception that, for 3=0, top surface strips Ab on both sides of the wing 
t^^tribute to the limit cycle. The following numerical data are presented 
for a typical oscillation cycle (cycle B in Figures 45 and 46): 
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1 


- 0.114 

U “ 301.1 m/sec 

dh » 0.00308 m 
5 ® 




Q * 18.2 KN/m" 
A » 0.0524 m 


» 4 ' 




-1 


^ 5= 0 . 3767 sec 

•7 

E * 2.82 m-N/cycle 

(of which 76.6% is due to mechanical damning. 21.4% is due to aerodynamic damping) 


Wk 




Ab AQ dh = 5.88 Ab m-N/cvcle 


PE » -I” “ 1.776 m-N/cycle 

Ab - 0.1425 X (~) « 0.077 m 


In evaluating the above data representing the cyclic and the random 

response cases simulating wing rock, the following conclusions may be drawn: 

1. Wing rock was generated by a limit cycle mechanism due to the fluct- 
uating pressure changes- on the wing top surface, especially near the 
wing tip region. In the cases observed, the major reason for the pres- 
sure changes was the motion induced Aa which alternately caused 
leading edge stall and recovery. The amount of energy input to the 
limit cycle, to be consumed by aerodynamic and/or mechanical damping, 
may be evaluated quantitatively. 

2. Based on kinematics, the roll oscillations acted together with a , 

o 

8 to cause substantial Aa changes which contributed to the limit cycle. 
Wing rock amplitude was thus dependent on the steady state sideslip 
angle 8 . 

3. With a sideslip angle 8, the upper wing surface on the windward side 
contributed the majority of the fluctuating pressures which caused 
wing rock. 
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With Aa,the local fluctiiatlng pressure coefficient followed the ramp curve 
during leading edge stall and the consequent loss of lift. Dynamic 
overshoot was observed in both the lift loss and recovery phases. 

5. Based on the observations of the wind tunnel test data, a preliminary 
mathematical model of aircraft wing rock was developed which illustrated 
the motion-coupled limit cycle mechanism. Since the limit cycle de- 
fining the energy input to the aircraft was dependent on a number of 
flight condition parameters, no uniqueness condition can be established 
at this time to define and predict wing rock. Nevertheless, it is known 
tnat wing rock was directly related to the lift loss (ramp) of the 

C vs. a curve covering the fluctuating- pressure- active region of the 
P 

wing. Additional acquisition and processing of controlled 
test data as those described above may further refine and quantify 
the preliminary mathematical model of aircraft wing rock described in 
this report . 
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Section 4 


SPECTRAL AND STATISTICAL PROCESSINGS OF TEST DATA 

In the previous section, the real-time pressure and response data as 
well as the static and dynamic pressure distributions were presented. Assuming 
that the random pressure and response date were stationary within a 
limited time span, standard spectral processing techniques may be applied. 

The spectral processing consisted of the determination of the power spectral 
densities of the pressure and response data, as well as the spatial correla- 
tion functions of any two sets of data. A brief outline of the processing 
techniques may be found in Reference 1. 

4.1 Typical Pressure Power Spectra - Phase I 

The dynamic pressure data acquired during Run 8, Phase I were processed 
to generate the power spectral density functions. The corresponding pressure 
distributions on the right wing upper surface were given in Figure 10. 

In Figures 47 and 48 , four sets of data were used to generate the power spectral 
density functions of pressure transducer Nos. 1 and 2, respectively. These corre- 
spond to a = 0 degree, a = 8 degrees, a = 12 degrees, a = 16 degrees. The 
data time span used in the processing were as follows: 


Freq. Range (Hz) 

Record Processed 
(Sec. ) 

No. of Records 
Used 

Total 
Processed 
Data Time 
( Sec . ) 

0 - 200 

2 

16 

32 

200 - 2,000 

.200 

16 

3.2 

2,000 - 10,000 

.040 

16 

0.64 


Referring to Figure 47, the PSD scale on the right side, and the fre- 
quency scale on the top of the plots were based on actual tunnel test 
conditions. The corresponding scales on the left and bottom sides of the plots 
were transposed scales of the PSD data if the tests were performed under 


I 
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an equivalent flight test condition of F-5A, The scaling conversion was 
based on the Q, V, and c data of F-5A test aircraft N6009, Flight No. 871 
(see Reference 1, pp. 43-55): 

M = .925 h = 10,668 m 

Q = 14.36 KN/m^ V = 274.3 m/sec. 

c = 2.356 m = 4“/12° 

= 17.24 X 10® 

In Figures 47 and 48 , for a = 0 degree and 8 degrees, the peaks at approx- 
imately 3000 Hz and above (tunnel test frequency scale) were identified as tunnel 
noises. (See Reference 5). These peaks became l^ess significant as the angle-of- 
attack and PSD level increased (a = 12°, 16°). The substantial increase in the 
PSD level at cy = 12 degrees signaled the local flow separation. 

For Of = 10 degrees, the pressure PSD data were processed for trans- 
ducer Nos. 2, 3, 4 using the same processed data time as shown in the 
table relating to transducer No. 1 pressure data. The results are 
presented in Figure ^9. At o? = 10 degrees, transducer Nos. 3 and 4 
were under the separated flow, while No. 2 was still in the unseparated flow 
region. The PSD level difference between the unseparated and separated flows 
was 10 dB or more except at the high frequency region. Because of the close 
proximity of the transducers, the similarity of the PSD patterns suggested a 
high degree of spatial correlation of the dynamic pressures. This point will 
be investigated further. 

In order to observe the influence of the wake flow to the horizontal 
tail during aircraft buffet, the pressure PSD data were processed for 
transducer Nos. 12, 22, 24 corresponding to o' = 0 degree, 8 degrees, 12 
degrees, 16 degrees. The results are presented In Figure 50. Transducer 
No. 12 was located at W.S. 71.4 or 47 percent semi-span and 40 percent 
chordwise location of the upper wing surface (see Figure 3), Transducer 
Nos. 22, 24 were located on the horizontal tail at the same semi— span loca- 
tion downstream of No. 12 and 40 percent of the local chord. No. 22 was on 
the upper surface and No, 24 on the lower surface. In examining the figure, 
it was noted that unseparated flow appeared at transducer No. 12 when 
Qi = 0 degree, 8 degrees, and also at transducer No. 22 when a = 0 degree. 
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Figure 30. Pressure PSD Plots of Transducer Nos. 12, 22, 24 for F-5A Scale Model, Rui 
Phase I, P= 0°, M = .926, <5^/ <5^ = 5°/12'^, 6^ = 0®. The Scales have been 
Modified to Conform to Flight Test 871 Conditions. 



Otherwise, the local flow was separated. The similarity of PSD curves of 
transducer No. 12 and No. 22 at ck' = 0 degree suggested a certain degree of 
spatial correlation, most possibly due to the dominance of characteristic 
wind tunnel noise. At cK = 8 degrees, 12 degrees, transducer No. 22 and No. 

24 featured almost identical PSD curves, reflecting a strong wake flow con- 
dition. At O' = 16 degrees, the very high PSD level of transducer No. 22 
indicated the high turbulence due to both wake flow and horizontal tail top 
surface flow separation. This PSD curve was one of the highest broad-band 
spectral level recorded in the test program. The high dynamic pressure 
intensity on the tail surface for large q^'s was expected to affect the control 
characteristics of the aircraft, as compared to the direct buffet pressure 
excitation through the aircraft wing. In the following subsection, the effect 
of flap settings will be described. 


4.2 Effect of the Flap Settings 

Effect of the flap settings was examined using test data of Run 8 and 
Run 39, Phase I. Run 39 was conducted in test model configuration 9. The 
configuration featured wing- tip-mounted Sidewinder missiles. The leading 
and trailing edge flaps were in the retracted positions = 0°/0°). 

The aileron and horizontal tail surface settings were 0 degree. The nominal 
.sideslip angle was 0 degree. The run was conducted at M = 0.925 with a dynamic 
pressure of 28.30-28.68 KN/m^ (591-599 psf ) . The nominal Reynolds number was 
11.483 X 10^/m (3.50 x 10^/ft). 
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A comparison of dynamic pressure distributions of Runs 8 and 39 at critical 

Ot’ s is shown in Figure 51. The major difference between Run 8 and Run 39 was 

the extended flaps of the former ( ~ 5°/12°) and the completely retracted 

flaps of the latter (the camber effect). Furthermore, there existed a difference 

f\ A 

in Reynolds number, 7.402 x 10 /m for Run 8 vs. 11.483 x 10 /m for Run 39. Since 
the two Reynolds numbers were of the same order of magnitude, the major difference 
in the pressure distribution may be attributed to the camber effect. (This 
point was borne out by detail comparison of the dynamic pressure data presented 
later.) ' 

In Run 39, as the angle-of-attack increased from 0 degree to 8 degrees, 
the separation-inducing shock moved speedily toward the leading edge in the 
wing-tip region. At a = 8 degrees, high rms dynamic pressure was observed 
at 85 percent semi-span close to the leading edge. With further a increase, 
the peak dynamic pressure location moved rapidly along the leading edge toward 
the fuselage. This behavior was different from the case of Run 8 where the 
shock-induced flow separation region appeared first throughout the semi-span 
in the trailing edge region at approximately a = 8 degrees. 

With additional increase in a, the shock front and the separated flow 
boundary moved upstream toward the leading edge. At a = 14 degrees, 16 degrees, 
the complete upper wing surface was subjected to the separated flow. The rms 
dynamic pressures were highest in the region from the fuselage up to 72.8 
percent semi-span, with the wing-tip region showing lower dynamic pressures. 

In general, at moderate angle-of-attack, the deflected leading edge flap 
of Run 8 tended to push the flow separation boundary further downstream as 
compared to the case of the retracted leading edge flap. The extended trail- 
ing edge flap tended to keep the shock-induced flow separated in region further 
inboard. This can be attributed to the increased camber in the inboard portion 
of the semi-span for the extended flap case. As far as the shock and the shock- 
induced flow separation region are concerned, the flow behavior at various wing 
sections was more uniform along the span-wise direction with the deflected 
flaps. The above-mentioned differences were most prominent at a = 6 degrees, 

8 degrees. With further increase in a, the differences in the flow behavior 
became less distinct as the upper surface was completely under separated flow. 


103 



5V12 





The pressure power spectral density plots of Run 39 are presented in 

Figures 52 to 54 for Transducers Nos. 2, 12, 18 respectively. Each set of PSD 
plots has four curves corresponding to = 0 degree, 8 degrees, 12 degrees, 

16 degrees. The scales have been modified to conform to Flight 825, Run 5 
condition of a F-5A test aircraft (see Reference 1). The modification of 
scales was consistent with the PSD data of Run 8 presented previously 
(Figures 47 to 50) . 

In comparing the PSD plots of transducer No. 2, Run 39 (Figure 52) with 
those of Run 8 (Figure 48), it was noted that the major difference in the 
PSD trend appeared at = 8 degrees. At this angle, the PSD level of Run 39 
at transducer No. 2 was substantially higher than the corresponding PSD of 
Run 8. This same peak may also be observed in the pressure rms plot in 
Figure 51 for jy =* 10 degrees. For the other qt angles (0®, 12°, 16°), the 
PSD plots of transducer No. 2 for the two runs were almost identical to 
minute details. In general, under otherwise identical test conditions, the 
difference in the flap angle settings seemed to affect the movement of the 
shock boundary, the peak dynamic pressure level where the shock appeared, 
and the resulting development of the separated flow. At a given station 
away from the flaps, the dynamic pressures were pretty much defined by 
whether the local flow was unseparated or separated and the instantaneous 
angle~of~attack, irrespective of the flap angle settings. 

In examining the PSD plots of transducer No. 12 which represented a 
typical mid-span location on the top surface of the wing (Figures 50, 53), 
it was noted that for the same the differences of PSD trends between Run 8 

and Run 39 were not substantial. At any frequency band, the difference in 
PSD levels for the two runs remained mostly at 5 dB or less. 

Figure 54 shows PSD plots of transducer No. 18 of Run 39 which was 
located on the top surface of the undeflected trailing edge flap. No 
comparable data were processed for Run 8. In general, the PSD' s of 
transducer No. 18 were similar to those of transducer No. 12 for the same c^*s. 
The exception was the PSD of transducer No. 18 at (^ = 16 degrees which 
resembled a white noise distribution throughout the frequency band of 
interest. 
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4. 3 Comparison of Pressure PSD Data for the 
Fixed and Flexibly-Supported Models 

Run 5, Phase II was conducted with the same model configuration (No, 2) 
as Run 8, Phase I, except that in the second test phase, the model was 
supported on a flexible system. The tip-mounted Sidewinder missiles were 
included in the configuration. The leading edge and trailing edge flap 
settings were 5®/12®. The horizontal tail surface and aileron settings 
were 0 degree. The nominal sideslip angle was 0 degree. The run was 
conducted at M = 0.925. The dynamic pressure was 18.2 KN/m^ (377.7 psf). 

The Reynolds number was 7.38 x lO^/m (2.25 x 106/ft ). The critical angle- 
of-attack was 10 degrees at which the model oscillated along the flexible 
roll axis. Typical dynamic pressure coefficient distributions on the right 
wing upper surface for Run 8, Phase I and Run 5, Phase II were given in 
Figure 22. The plots cover the angle-of-attack range = 6 degrees, 10 
degrees, 12 degrees. 

In Figure 22, at ^ = 6 degrees, the rms pressure coefficient distribu- 
tions for Run 8, Phase I and Run 5, Phase II are qualitatively similar. Both 
reflected the development of the shock- induced flow separation covering approxi- 
mately 30 percent of the chord toward the wing trailing edge. At jy = 10 
degrees, when roll oscillation took place for the flexibly-supported model, the 
Phase II results indicated a somewhat less distinct shock front location, 
possible due to the oscillatory change in the effective angle-of-attack caused 
by the wing roll motion. For = 12 degrees, the complete wing upper surface 
was under separated flow. The oscillating model featured more evenly distributed 
rms pressure coefficients as compared to the Phase I data. 

Examination of the static pressure distributions for Run 8, 

Phase I and Run 5, Phase II (not shown in this report) showed that 
the Cp coefficient developments on the wing upper surface with a - changes 
remained essentially the same for the two runs. In either case, the highest 
Cp amplitude (negative pressure) developed at approximately a = 10 degrees, 
when the flexibly-supported model started to oscillate in roll. The C 
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amplitudes showed moderate reduction as the angle-of-attack was further 

increased. The upper surface of the horizontal tail exhibited moderate 

negative (in order of —0.10) up to = 12 degrees. Beyond which, the 

amplitude increased substantially with At ^ = 16 degrees, the C at 

P 

transducer No. 22 reached —0.7142 in Run 8, Phase I and —0.7062 in Run 5, 

Phase II. 

In Figures 55,56, the power spectral density functions of pressure 
transducer Nos. 1, 12, 22, 24 for Run 5, Phase II are plotted with those 
corresponding to Run 8, Phase I shown previously in Figures 47,50 for comparison 
purpose. The angle-of-attack range covers a = 0 degree, 8 degrees, 12 degrees, 

16 degrees. The scales used were modified to conform to the equivalent flight 
test condition of the F-5A, as was the case in Figure 47. In both figures (55,56), 

two arrowheads were attached to the PSD plot at the upper left corner. They 

identified the pressure PSD peaks corresponding to (1) the sting bending for 
the supporting systems, with a natural frequency of approximately 6.6 Hz, or 

0.865 Hz after the frequency scale was modified to conform to Flight Test 871 

conditions, (2) the first wing symmetrical bending, with a natural frequency 
of approximately 53.5 Hz, or 7.01 Hz based on the modified frequency scale. 

Referring to Figure 55 , the PSD data for transducer No. 1 and No. 12 
were almost identical for the two runs with the exception of the very low 
frequency region corresponding to the roll oscillation natural frequency of 
the flexible support. The slightly higher PSD level of the flexibly- 
supported model gave a good indication on the effect of the induced dynamic 
pressures due to model oscillation. The peak at approximately 340,6 Hz (based 
on the modified frequency scale) or 2600 Hz in actual test condition was 
characteristic PSD peak of the tunnel noise. The tunnel noise peak was pro- 
minent at low angle-of-attack when the overall dynamic pressure levels on the 
wing surface were low. 

The corresponding tail surface pressure PSD comparison for the two runs 
is shown in Figure 56 for transducer No. 22 and No. 24. In this case, more 
prominent deviations of the PSD levels were observed, especially at = 8 
degrees. It is reasonable to assume that the differences at q/ = 8 degrees, 
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Figure 56. Comparison of Horizontal Tail Surface Pressure Data for Configuration 2 for the Regularly- 

Supported Model (Run 8, Phase I, Indicated in Dotted Lines) and the Flexibly-Supported Model 
(Run 5, Phase II, Indicated in Solid Lines), P = 0°, M = 0.925, <5^/6^ = 5°/12°. 6^ = 0®. 

The Scales are Modified to Conform to Flight Test No. 871 Conditions. The Arrowheads Identify 
(A) the Sting Bending Mode and (B) the First Wing Symmetrical Bending Mode. 








prior to large amplitude model oscillation (a* = 10°), were due to minute test 
condition changes which affected the dynamic pressures in the wake flow region. 
At = 16 degrees, the PSD curves for the two runs were almost identical, as 
were the case of the wing surface transducers (Figure 56). 


4.4 Comparison of Measured Power Spectral Data 
with Flight Test Spectral Data - Reynolds 
Number Effect 

Comparison of measured power spectral data with equivalent flight test 
spectral data for M = 0.925 case was made using the data recorded during Run 
39, Phase I test. The features as well as the development of flow pattern 
of Run 39 were described in subsection 4.2. The right wing top surface static 
and dynamic pressure coefficient distributions are plotted in Figure 57 . 

Figure 58 are the pressure PSD plots for transducer No. 2 of the scale model 
during Run 39, Phase I test. In these plots, only converted PSD and frequency 
scales (to Flight 825 conditions) were presented. The three plots consisted 
of PSD data as described below: 


a = 8° 

marked 

(a) 

a = 12° 

marked 

(b) 

a = 16° 

marked 

(c) 


All scale model PSD data were plotted in solid lines. For comparison 
purposes, the corresponding PSD data of a transducer at the identical geo- 
metrical position on the F-5A test aircraft acquired in Flight 825 Run 5 

(M = 0.925, h = 10,668m, 6 /6^ = 0°/0% K, = 17.24 x 10^) were overplotted 
o n t In 

on the figure in broken curves. Specifically, the PSD data were based on 
approximately 0.5 second pressure records of a transonic maneuver performed 
by the test aircraft. These PSD data were identified as: 

a =7.2° -8.6° marked (A) 

a = 11.4° - 13.2° marked (B) 

CL = 13.9° marked (C) 

No comparison was made at the beginning of the maneuver. 
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Figure 58. Pressure PSD Plots of Transducer No, 2 for Scale Model in Run 39, 

Phase I and F-5A Aircraft in Flight Test 825, Run 5, M = 0.925, 

“ 0°/0°, 6, = 0°. The Scales are Based on Flight Test 825 Conditions 





In comparing the scale model and the aircraft pressure PSD data at trans- 
ducer location No. 2, it was noted that the PSD levels and the frequency dis- 
tributions of the tunnel test and flight test data were quite comparable. 

For a typical inboard transducer. No. 12 of Figure 3, the pressure PSD 
data are presented in Figure 59. The identification letters for the scale 
model and the F-5A aircraft flight test data conventions were the same as in 
Figure 58, In comparing the scale model and the aircraft pressure data, it 
was found that except at a = 8°, the PSD levels and the frequency distributions 
of the tunnel test and flight test data were again quite comparable. At 
a = 8®, a difference of PSD level as much as 15 dB was observed. The peak 
of the tunnel test data at approximately 230 Hz and above were caused by 
tunnel noises. For the scale model, the PSD levels also showed an average 
of 20 dB increase when a was changed from 8® to 12°, indicating a drastic 
expansion of the separated flow region on the wing surface in the inboard 
area. This phenomenon was also observed in Figure 53. In comparing the 
two cases for station 12, it seemed that in the transition stage (from sep- 
arated flow to shock-induced flow separation) , the actual aircraft local 
transition appeared at a lesser angle-of-attack. After the transition, the 
local pressure rms level corresponding to a separated flow condition were 
about the same for the actual aircraft and for the scale model. 

Similar comparison was made for the flight and model test results for 
M = .75 case. The test data of Run 43, Phase I were used. Run 43 was con- 
ducted with model configuration 9, i.e., all flaps and control surfaces at 

retracted positions. Two sidewinders were mounted on the wing-tips. Run 43 

2 

was executed at M = 0.75. The dynamic pressure was 30.23 KN/m (631.3 psf ) . 

6 6 

The Reynolds number was 13.99 x 10 /m (4.264 x 10 /ft.). The model Reynolds 
number was 4.71 x 10^. At M = 0.75, the flow separation was leading edge 
induced. The phenomenon may be traced through the static and dynamic pressure 
coefficient plots with increasing a (Figure 60). At a = 6 degrees, separation 
existed in 85 percent semi-span location up to the tip. At a = 10 degrees, 
the separation region moved inboard toward 33 percent semi-span location. 
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At (y = 12 degrees, the complete upper wing surface was subjected to separated 
flow. This separation behavior was similar to the F— 5A aircraft test data that 
also showed gradual development of flow separation with increasing at 
M = 0.75 (Figure 34, Reference 1). 

The PSD plots of pressure transducer Nos. 2, 5 and 11 based on Run 43 
and a “ 8 degrees, 12 degrees and 14 degrees are presented in Figures 6] to 63, 
respectively. The Run 43 data are identified by lower case letters a, b, and 
c. The PSD and frequency scales of the scale model data are on the righthand 
side and the bottommost portion of each plot. Also plotted on these figures 
are the corresponding PSD data obtained in Flight Test 825, Run 7 under the 
following flight conditions: 


M = 0.75 

h = 7,772 m (25,000 ft, ) 

6n/6f = 0V0° 

= 18.90 X 10^ 

Store: Wing- tip-mounted Sidewinders 

Run 7, Flight 825 aircraft configuration and the Mach number were identi- 
cal to the corresponding scale model test conditions of Run 43. Run 43 had 
a zero sideslip angle, while the sideslip angle of Run 7, Flight 825 was 
variable throughout the maneuver. The flight test PSD data are identified by 
capital letters A, B, and C, with the appropriate angle-of-attack range marked 

r 

under each plot. The flight test and scale model data are plotted over each 
other taking into consideration the model scaling relation. The flight test 
PSD and frequency scales are on the left side and at the bottom directly 
under each plot. The arrowheads on the pressure PSD plots (Figures 6l to 63) 
for = 8 degrees identify the peak frequencies corresponding to the rigid 
body and flexible mode motions of the scale model as defined in subsection 
2.2 (page 13) and Figure 4 (page 17). 

An examination of Figures 61 to 63 indicates that with allowance made for 
the difference in the flight test and model test conditions including the 
sideslip angle, the continuous changing angle-of-attack during flight, the 
difference in the structural modes of the aircraft and model, etc., the 
developments of the dynamic pressures for the two cases could be correlated. 
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Figure 62. Pressure PSD Plots of Transducer No. 5 for Scale Model in Run 43, Phase I and F-5A 
Aircraft in Flight Test 825, Run 7, M = .75, 6 /d = 0°/0®, 6 = 0°. The Left- 
Hand Scale and the Upper Frequency Scale are Based on the Flight Test Data. The 
Right-Hand Scale and the Lower Frequency Scale 'are Based on the Scale Model Data. 
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Figure 63. Pressure PSD Plots of Transducer No. 11 for Scale Model in Run 43, Phase I and 
F-5A Aircraft in Flight Test 825, Run 7, M = .75, 6 /6 . = 0^/0% 6 = 0°. The 
Left-Hand Scale and the Upper Frequency Scale are Based on the Flight Test Data. 

The Right-Hand Scale and the Lower Frequency Scale are Based on the Scale Model Data. 




In general, the flight test pressure PSD levels were somewhat lower (in 
the order of 5 dB) than the corresponding model test data. The consistent 
deviation in the spectral level might be contributed, at least partially, 
to the Reynolds number effect corresponding to the flight and tunnel test 
conditions (R^ = 18.90 x 10^ vs. 4.71 x 10^), 

In the transonic range with flow separation, it is known that the Rey- 
nolds number effect is significant for wings with thick airfoil, such as those 
with a thickness/chord ratio of 8 to 10 percent. In actual flight conditions 
where the Reynolds number is high, the shock-induced transition point is further 
downstream as compared to the low Reynolds number case if other conditions are 
identical. In model tests for the same flow condition, except with a lower 
Reynolds number, the flow transition strip tends to trigger the flow separation 
prematurely. The result is less favorable lift distribution and drag. Further- 
more, the location of the transition strip is known to affect the flow separation 
significantly for thick airfoils. 

For thinner airfoils - e.g., the F-5A with NACA 65A004.8 wing airfoil 
sections - it is expected that the Reynolds number effect Is less significant. For 
M = 0.75 and 0.925, the range of Reynolds numbers of the available data is as 
follows : 


TYPE OF TEST 

M = 0.75 

M = 0.925 

Flight Tests 

18.90 X 10^ 

17.24 X 10^ 

Tunnel Tests 

2.48 X 10^ to 4.71 X 10^ 

2.48 X 10^ to 4.71 X 10^ 


The flight test and model test data for the F-5A were examined to evaluate 
the Reynolds number effect. Since only dynamic pressure transducers were used 
in the flight tests, which did not yield valid static pressure data, a 
comparison of the dynamic pressure coefficients and spectral makeups in the 
separated flow region as well as the transition line as a function of a were 
used to determine the Reynolds number effect. 

Figures 64 and 65 present the dynamic pressure coefficient distributions 
of Run Nos. 39 and 43, Phase I test and corresponding flight test 825 data 
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for a = 4®, 6°, and 8°. The model Reynolds numbers for Run Nos. 39 and 43 

6 6 

were 3.87 x 10 and 4.71 x 10 , respectively. 

Based on the data shown in Figures 64 and 65, it was found that the 
trend of development and expansion of the flow separation region were similar 
for the wind tunnel model test and the corresponding flight test. In general, 
measured in angle-of-attack, it took the scale model approximately an addition 
2®-4° to develop a similar separation flow pattern on the top wing surfaces. 

For M = 0.925 case (Figure 64), at a given spanwise location, the shock tended 
to stay closer to the trailing edge on the scale model as against the flight 
test results. This observation appears to be contrary to previous test 
results. For instance, References 2,3 show that under identical conditions 
except for Reynolds number, the shock on the scale model always appeared 
further upstream in relation to the full-scale aircraft wing section when the 
airfoil section thickness ratio was high. The phenomenon observed in the pres- 
ent test was somewhat perplexing since transition strips were installed on this 
F-5A scale model which should trigger the shock-induced flow separation in 
an upstream location on the wing chord; i.e., the artificial boundary layer 
transition together with testing at less than flight Reynolds number should 
result in a relatively thicker model boundary layer, hence, leading to earlier 
shock- induced separation. In evaluating this problem, certain factors 
were noted which were believed to contribute to the scale model flow behavior; 
(1) the aeroelastic effect of the full-scale aircraft caused an increase in 
the local angle-of-attack in the outer span of the wing. It was known that 

in F-5A aircraft, under the conditions where the pressure data were evaluated 

2 

(M = .925, h = 10,668m, Q = 14.36 KN/m ), the wing-tip was twisted in such a 
manner that its local angle-of-attack was increased ly to 2°; (2) the transient 
effect in a transonic maneuver of the aircraft where the angle-of-attack 
changed continuously; (3) F-5A wing has a thin airfoil section, namely 
NACA 65A004.8. 
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4.5 Typical Response Power Spectra and Dynamic Aeroelastic Effect 

It was well known that under moderate to high dynamic pressures, a fairly 
rigid scale model such as the F-5A model will respond to the tunnel aero- 
dynamic forces with measurable displacements. With this point in view, the 
F-5A scale model was instrumented to measure certain accelerations and wing 
root section moments as described in Table 1. The dynamic responses were 
dependent on the model vibration modes and the structural characteristics of 
the supporting system, The modal data of the fixed and the flexibly- supported 
models were described in subsection 2.2. Even though these dynamic charac- 
teristics were quite different from the real aircraft, the model response 
data may be analyzed and extended to predict full-scale aircraft responses 
under dynamic buffet pressures. 

For Run 8, Phase I the right wing-tip acceleration (la) data were pro- 
cessed into power spectra corresponding to q/ = 0 degjree, 8 degrees, 12 degrees, 
16 degrees. The spectral processing was limited to 0—200 Hz using the average 
of 16 records, each of 2 seconds duration. The processing bandwidth was 1/2 
Hz. The power spectra obtained in this manner were presented in Figure 66. 

The PSD and frequency scales were the test data scales. For the power spectrum 
corresponding to = 0 degree, certain peaks were identified with letters 
described as follows (see Figure 4): 

A. Sting and Balance Bending, 6.6 Hz 

B. Balance Roll, 13,6 Hz 

C. Model Wing 1st S)mimetrical Bending, 53.5 Hz 

D. Rigid Roll Plus Wing Anti- Symmetrical Bending, 71 Hz 

E. Wing Anti- Symmetrical Bending, 82,7 Hz 

F. Sting and Support System Torsion, 105 Hz, The peaks appear at 
120 Hz in some PSD plots may be due to the harmonics of the 
line noise, 

G. Wing 2nd Symmetrical Bending, 171 Hz. 
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In examining Figure 66 , it was noted that in the very low frequency 
region (under 10 Hz), the increase per degree of angle-of-attack in response 
spectra (and rms level) was moderate for Ot = 0 degree to Ot = 8 degrees, when 
compared with corresponding increase for a = 8 degrees to a = 12 degrees. In 
reviewing the static and dynamic pressure data (Figure 10) , a = 8 
degrees may be considered as the approximate condition where buffet onset 
took place. The model responded with increasing rate of intensity after 
buffet onset. It was also noted that except for the mode corresponding to model 
support system bending, the change in wing- tip responses from ^ = 12 degrees 
to (y = 16 degrees was very moderate. This was consistent with the dynamic 
buffet pressure data, where the spectral density level reached a plateau cor- 
responding to large angles- of —attack well after buffet onset (See Figures 
47., 48). 

In Figure 67 ,b 8 , the acceleration PSD's at the left wing— tip (2a) and the 
CG (3a) were processed and plotted. The test conditions were identical to 
those of Figure 66 . The peaks of the PSD's at (y = 0 degree were identified 
with the same letters as in Figure 66 . A comparison of Figures 66,67, showed 
that the left and right wing-tip responded to both symmetrical and antisym— 
metrical modes with fairly even intensity for all angles-of-attack. The only 
exception was for peaks (A), (B) at o' = 16 degrees. The difference in the tip 
responses for the latter case seemed to indicate some asymmetrical model motion, 
possibly due to the asymmetrical mass inertia distribution (cut-outs, channels, 
etc.) in the wing structures. The CG PSD plots (Figure 68 ) showed the ever 
increasing PSD intensity with ^ corresponding to the support system bending (A). 
As was expected, the peaks corresponding to the model antis 3 nranetrical modes 
(D), (E ) were relatively low at the CG. 

The roll moment data acquired by the 6 — component balance were processed into 
PSD' s and presented in Figure 69 . Peak identifications for the PSD of 
= 0 degree were the same as the previous figures (Figures 66 *“ 68 ) , with 
the exception that the peak at (C) (60 Hz) seemed to indicate a line noise. The 
120 Hz peak at (P) also contributed to line noise harmonics. 
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Figure 69. F-5A Scale Model Roll Moment PSD Plots, Run 8, Phase 
M = .926, 6 = 5V12\ 6, = 0 ° . 






For the F-5A scale model, semi-conductor type strain gages (Kulite 
ULP-120-160) were installed near the semi-span root sections. Since the wing 
structure was made of solid steel with cut-outs and harness channels, etc, , it 
was not possible to install gages which yielded sectional bending and torsion 
moment data directly. Through a pretest calibration procedure using various 
loading conditions, calibration formulas were developed to give approximately 
sectional bending and torsional moments based on gage outputs. Applying 
these formulas, the bending and torsional moments were acquired for Run 8. 

The corresponding PSD plots for four angles-of— attack (q/ = 0®, 8® 10®, 12°, 

16®) were presented in Figures 70, 71, Identical peak mode designations (as 
in Figure 66) were applied on these figures for (y = 0 degree. Referring to 
Figure 70, the bending moment PSD level behavior as a function of (y was some- 
what similar to the wing-tip acceleration PSD level variations (Figure 66,67). 

The torsional moment plots (Figure 71) featured lower PSD levels as compared 
to the bending moments. This was to be expected since the F— 5A wing featured 
a moderate sweep back. Again the trend of PSD level increases with a was consis- 
tent with the dynamic pressure and acceleration PSD data presented previously. 

In order to demonstrate the impact of the rigid body and flexible mode 
motions to the unsteady pressure forces, the model right wing-tip acceleration 
PSD for Run 43, Phase I, at (^ = 8 degrees is presented in Figure 72. Refer- 
ring to the figure, the spectral peaks noted as (A) through (G) were 
identified previously. The seven peak frequencies were identified by arrow- 
heads on the corresponding pressure PSD plots of Figures 61-63. An examina- 
tion of the figures indicated that only corresponding to mode (A), the low 
frequency sting and balance bending mode, there appeared prominent pressure 
peak. The other modes (B) through (G) were not prominently represented in the 
corresponding pressure PSD plots. In general, it might be concluded that 
only very low frequency modes (10 Hz or less) had a prominent effect on the 
separated flow dynamic pressures. 
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Figure 70. F-5A Scale Model, Right Wing Root Section Bending Moment PSD Plot 
Run 8, Phase I, P = 0° , M = .926, 6 /d^ = 5®/12°, 6 = 0®. 
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Figure 71, F-5A Scale Model, Right Wing Root Section Torsional Moment 
PSD Plots, Run 8, Phase I, P= 0®, M = .926, 6 / 6 . = 5®/12 
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4.6 Low Frequency Range Pressure and Response PSD Data 
Corresponding to Maximum Wing- Rock Motion 

As described previously in subsection 3.6, during Phase II test, the 
model oscillated in roll as the angle-of-attack setting approached a critical 
value. In order to define the correlation between the roll oscillations and 
the dynamic buffet pressures or the model accelerations, PSD plots of fluc- 
tuating pressures and responses were generated to cover the lov? frequency 
range (0.5 Hz to approximately 100 Hz) for two typical cases when the max- 
imum wing-rock type motion occurred. For Run 26, a = a* = 10°, the model 
oscillated steadily and continuously throughout the recording period (Fig- 
ure 35). For Run 5, Ot = a* = 10°, the motion was highly nonstationary. 

The oscillations started in a random manner, reached a maximum amplitude and 
then died out (Figures 36) . Corresponding Phase I data (Run 51 and Run 8) 
were also processed for comparison purposes. The results are presented below. 

Run 26 of Phase II test was conducted in test model Configuration 11, 

6 /6 = 0°/0°, 6, = -10°, 6 = 0°/0°, 3=8°, Sidewinder missiles at wing- 

n f h a 

tips, M = .925, and Oi* = 10°. The corresponding Phase I test with a reg- 
ular sting mount was Run 51. Typical pressure and response data of both 
runs were processed into power spectra at Oi = 10° where the steady, high- 
amplitude roll oscillations occurred during Phase II test. The spectral 
processing was limited to 0.5 - 100 Hz using the average of 8 records, each 
of 4 seconds duration. The acceleration spectral data at the right wing- 
tip and the CG are presented in Figure 73. The plots on the left hand side 
are for Run 51, Phase I. The right hand side plots are for Run 26, Phase II. 
Referring to the response data, the major peaks were identified as below:- 

PHASE I PHASE II 

A Balance-Sting Bending Sting Bending 

B Balance Roll Roll Device Oscillation 

C First Wing Symmetrical Bending First Wing Symmetrical Bending 

The PSD plots for right wing-tip acceleration indicated a difference of 37 
dB in the first response peak (under 10 Hz) between the two runs, reflecting 
the large-amplitude roll oscillations in Phase II. 
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Figure 73. PSD of Right Wing-Tip Acceleration and CG Acceleration of the F-5A 
Scale Model, Run 51, Phase I (left) vs. Run 26, Phase II, a = 10°, 
B = 8°, M = 0.925, 6 = 0°/0°, 6^ = -10°. 
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In examining the pressure data» the dynamic pressure coefficient 
plots for Run 51, Phase I and Run 26, Phase II at a = 0° and a = 10® are 
shown in Figure 74. It was observed that the dynamic pressure levels in 

the inboard portion of the right wing upper surface for Phase II data at 
O' = 10® were much higher in comparison with the Phase I data. This seemed 
to be a common trend for the windward side of the wing semi-span of the 
model with a sideslip when a approached the critical angle-of-attack. 

The pressure PSD plots for transducer Nos. 3, 10, 11, 17, 22, 24 are presented 
in Figures 75 to 77 respectively. The locations of the transducers are given below: 


Transducer 

Number 

Surface 

% Semi-Span 

% Chord 

3 

RH Wing Upper 

85 

40 

10 

RH Wing Upper 

61 

21 

11 

RH Wing Upper 

61 

60 

17 

RH Wing Upper 

33 

42 

22 

RH Tail Upper 

47 

40 

24 

RH Tail Lower 

47 

40 


For Run 51, Phase I, the pressure peak corresponding to the Balance-Sting 
bending mode (as identified by letter A in the left hand side plots of 
Figures 75 to 77) may be clearly observed, while the pressure peaks corres- 
ponding to the other modal response frequencies above 10 Hz, if observable, 
were obscure. On the other hand, the pressure PSD' s for Run 26, Phase II, 
featured a prominent peak corresponding to the natural roll frequency (as 
identified by letter B in the right hand side plots of Figures 75 -77), with 
a number of unidentified lesser pressure peaks below 50 Hz, Corres- 
ponding to the first wing bending mode (51.6 Hz), no prominent pressure 
peak was in evidence. It was noteworthy that the pressure peak level at 4 Hz 
was extremely high for most of the transducers. No corresponding flight test 
results have yielded such a pressure level during wing rock in transonic flight. 

A key reason for the lack of a corresponding high pressure peak for aircraft 
wing rock might be due to the lack of sustained low frequency oscillation data 
covering a long time span. The pressure spectral data at frequencies above 
100 Hz, not presented here, showed much less differences in both amplitude and 
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F ■ gure 76. Comparison of Wing Pressure PSD Data for tlie Fixed Model 
Run 51, Phase I (left) and the Osc il lating M odel . Run 26 
Phase II, a = 10°, B = 8°, M = 0.925, 6^/6 j- = OVO", = 


I 
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Figure 77. Comparison of Horizontal Tail Surface Pressure PSD Data for the 
Fixed Model, Run 51, Phase I (left) and the Oscillating Model, 
Run 26, Phase II. a= 10°. 3 = 8°, M = 0.925,6 /6. = 0°/0°, 6. 








the spectral make-up for the two runs except for the transducers located on 
the inboard portion of the wing semi-span. The data presented above demonstrated 
the frequency-dependency, as well as the dependency on the deflection modes, of 
the dynamic pressures as excited by the flexible responding structures. 

For Run 5, Phase II test, the random appearance of the roll oscillations 
and the subsequent reduction in oscillation amplitude made the recorded 
pressure and response data highly nonstationary. In contrast, in order to 
maintain a high degree of frequency resolution as required in the low fre- 
quency range, and to keep the standard error to an acceptable level, sub- 
stantial duration of stationary real time data are needed. It is well 
known that the normalized standard error e for a highly resolved PSD 

A 

estimate G (f) from a sample record x(t) with a true PSD function G (f) is 

X X 


s,d, G (f) 


(B tT 
e 


where = the equivalent bandwidth of the filter 
T = the sample record length 


By examining the real time data of the roll oscillations as shown in 

Figure 36, it was noted that the high amplitude oscillations lasted onlv 
2—4 seconds. The equivalent bandwidth of the filter (analog) used in the 
spectral processing was 0.375 Hz. In order to maintain reasonable statis- 
tical accuracy, Run 5, Phase I data were processed using three different 
sample record lengths (8 sec, 16 sec, and 32 sec) to cover the time span 
where roll oscillations were at high amplitudes. In Figure 78, selected 
spectral data of the model response (Right Wing-Tip Acceleration) and the 
unsteady aerodynamic pressure (Tranducer No. 3) are presented for comparison 
purpose. Plots (a) of Figure 78 are for Run 8, Phase I. Plots (b), (c) , 

(d) are for Run 5, Phase II. The sample record lengths used and the corres- 
ponding normalized standard errors are described below: 


(a) 

T = 32 

sec 

e = 

.2886 

(b) 

T = 32 

sec 

£ = 

. 2886 

(c) 

T = 16 

sec 

e - 

.4082 

(d) 

T = 8 

sec 

e = 

.5774 
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Figure 78. PSDs of Pressure No* 3 and Right Wing-Tip Acceleration Based on Run 8, 

Phase I (left) and Run 5, Phase II Data, a= 10", B = 0", M = 0.925, 6/6. = ! 










The major peaks of the response data were again identified by letters 
A, B, and C as in the previous case. For the pressure spectral data, only 
the peaks corresponding to the Sting-Balance Bending mode (A) for Phase I 
and the Natural Roll Frequency (B) for Phase II are identified. In exam- 
ining the Phase II plots, it was observed that the shorter the duration of 
the sample record length, the higher the peaks values. This was to be ex- 
pected. For the same reason, no prominent peak corresponding to the Nat- 
ural Roll Frequency was observed on the pressure PSD plots. 

Response data and pressure data at key locations recorded during Run 
5 and Run 26 of Phase II test and Run 8 of Phase I test at cx = 10° were 
digitized to determine the input-output phase relationship as well as the 
picture of the makeup of the rigid body modes (roll, yaw, pitch) and the 
mode coupling. The recorded data were digitized with a sampling rate of 521 
samples/second, and a low-pass analog filter with a cutoff frequency of 220 
Hz was used during the digitizing process. The data being digitized are 
listed below. 

Run 5, Phase II, a = 10° — Roll, yaw, right and left wing-tip 

accelerations (la, 2a), transducer 
Nos. 1, 3, 4, 12 

Run 8, Phase I, a = 10° — Right and left wing-tip accelerations 

(la, 2a), transducer Nos. 1, 3, 4, 12 

Run 26, Phase II, CL = 10° — Roll, yaw, pitch, right and left wing- 

tip accelerations (la, 2a), transducer 
No. 4 

The Fourier transforms of the above-listed data were generated by a 
Spectral Analysis computer program using the Fast Fourier Transform method. 

Four records, each of 2 seconds duration, were used for averaging purposes, 

(The segment used in Run 5, Phase II data is indicated by letter B in Figure 
36a. ) The Fourier transform data corresponding to two specific frequencies were 
of interest, i.e., the natural frequencies corresponding to the roll device 
oscillation for Phase II test (4.07 Hz) and the first wing symmetrical bending 
mode (51.88 Hz). The results on the relative phase relationship (in degrees) 
are presented below: 
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Run 5, Phase II, 
a =10°, M = 0.925 
= 5°/12°, 

■5,, = 0“, S = 0° 


Run 8, Phase I, 

OL =10°, M = 0.925 
<Sn/<Sf = 5°/12°, 

= 0 °, 6 = 


Run 26, Phase II, 
a =10°, M = 0.925 
(S /6 = 0°/0°, 

nr 


5 = - 10 ' 

h 


,3 = 8' 


Signal f=A.07 Hz f=51.88H: 
Ident. 9^-9 . 9-9 


f=51.88 Hz 


f=A.07 Hz f=51.88 Hz 


Pitch 


161.8' 


-14.4' 


177.0' 
-104.9' 
-98.8' 
161 .0' 
177.9' 


-30.7° 

-113.8° 

-118.7° 

-156.9° 

-109.7° 


-130.6° 

- 10 . 0 ° 

-33.3° 

-38.6° 

-37.6° 


170.4' 


The last row gives the phase angle'- between the modal responses (la) and 
the modal forces contributed by the pressures measured at the transducers 
listed above, corresponding to the first wing symmetrical bending mode. Based 
on the data obtained, it was clear that the yaw angle and pitch angle os- 
cillations had the same major frequency with a slight phase shift from the 
roll oscillations. 


4.7 Wake Flow Effect and Horizontal Tail Buffet 

Even prior to flow separation on the aircraft wing, the wake flow and its 
effect on the tail surfaces and the overall responses of the aircraft were very 
complicated. Wing surface flow separation (shock or otherwise induced) further 
complicated the picture. In addition, the tail surfaces might also encounter 
shock- induced flow separation, depending on the aircraft angle-of-attack and 
the tail surface setting. 

In order to examine the wake flow effect on the horizontal tail due to 
flow separation, the dynamic pressure data of Run 19, Phase I were processed. 

Run 19 was conducted in test model configuration 3. The configuration featured 
wing— tip-mounted Sidewinder missiles. The leading edge and trailing edge flap 
angles were 5®/12“. The aileron and horizontal tail surface settings were 
0 degree and -10 degrees, respectively. The nominal sideslip angle was 
0 degree. The run was conducted at M = 0.925. 

Run 19 covered the following angle-of-attack settings; 0°, 4®, 6®, 8®, 

10®, 12®, and 14®. Among these, the dynamic pressure data of transducer Nos. 
14, 17, 22, 23, 24, and 25, corresponding to = 0®, 8®, 10®, and 14° were 

digitized in order to perform spatial correlation for each pair of transducers 
using a spectral analysis computer program. Transducer No. 14 was located at 
47 percent semi— span and 80 percent chordwise location on the trailing edge 
flap upper surface of the right wing. Transducer No. 22 was located on the 
upper surface of the horizontal tail at the same semi-span location downstream 
of No. 14 and at 40 percent chordwise location. Transducer 24 was located on 
the lower surface, directly under transducer 22. Transducer Nos. 17, 23, and 
25 were located at 33 percent wing semi-span, with No. 17 at 42 percent chordwis 
location of the right wing top surface. Transducers 23 and 25 were at 40 per- 
cent chord of the horizontal tail with No. 23 on the upper surface and No. 25 
on the lower surface. The severity of the wake flow and its effect on 
tail buffet may be estimated through the following; 

1. Whether or not flow separation takes place in that part of the 
wing span directly upstream of the tail surface. 

2. High PSD level of the dynamic pressures on the tail. 
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3. Good correlation existed between dynamic pressures measured on the 

wing surface near the trailing edge and on the tail surface downstream. 



The dynamic pressure data measured at the above-mentioned transducers were 
digitized with a rate of 10,400 samples /second and the corresponding folding 
frequency was 5200 Hz. A low pass analog filter with a cutoff frequency of 
4500 Hz was used during the digitizing process. Spatial correlation between 
transducer pairs were performed in the form of cross-correlation, cross-power 
spectra, and coherency function. In order to show the general trend of the 
coherency function, a time span of 0,2 second was chosen. The maximum time 
lag used in cross-correlation computation was one fourth of the time span. 

Prior to presenting the spatial correlation data of the transducer pairs, 
the static pressure coefficient the rms dynamic pressure coefficient 

(C-) distributions on the right wing upper surface were examined to trace the 
shock /bound ary layer interaction and the separated flow regions. The results 
are presented in Figure 79 . It can be seen that transducer No. 14 was in the 
separated flow region for o' > 8°, while the local flow at transducer No. 17 
became separated at approximately a = 12® 

Spatial correlation calculations were performed for all possible pair 
combinations of pressures at the six transducer locations. High level correla- 
tions existed between (24, 25) for a = 0®; (14,22) and (24,25) for a = 10®; 

(22, 23) and (24, 25) for a = 14®, Typical spatial correlation plots for 
transducer pair (24, 25) at o' = 0®, (14, 22) at a = 10°, and (22, 23) at a = 14® 
are presented in Figures 80-82 > respectively. On plot (b) of those figures, 

the cross-power spectrum and the coherency function are presented. The unit or 

2 2 

the modulus is (N/m ) /Hz, the phase angle unit is degree. On plot (a), the 

auto-correlation functions and the normalized cross-covariance function (the 

cross— correlation coefficient) are given. The unit of the auto-correlation 

2 2 

functions (the top two plots) is (N/m ) . In addition to spatial correlation, 
cross-power spectra, etc,, the power spectral density functions of the six 
pressure transducers were also processed based on the digitized pressure data 
corresponding to various ang le-of-attack settings. For each power spectrum, 
data were processed using three continuous frequency ranges with different band- 
widths, The amounts of real time data used for the three frequency ranges are 
described in the following table. 
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(a) The Auto and Cross-Correlation Functions 

Figure 80. Spatial Correlation for Transducer Pair (24, 25) of 

Run 19, Phase I, rv = 0°, P = 0°, M = 0.925, 6^/6^ = 5V12°, 
6^ = -10° (Continued) 
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(b) The Cross-Power Spectrum and the Coherency Function 

FlEure 81. Spatial Correlation for Transducer Pair (14, 22) of Run 19 
Phase I, « = 10", P = 0", M = 0.925, 6^6^ = 5°/12°, 

6, = -10° (Concluded) 
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(a) The Auto and Cross-Correlation Functions 

Figure 82. Spatial Correlation for Transducer Pair (22, 23) of Run 19, 

Phase I, „,= 14“, p = 0”, M = 0.925, 6^/6^ = 5V12 , 6^ = -10 

(Continued) 
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(b) The Cross-Power Spectrum and the Coherency Function 

Figure 82. Spatial Correlation for Transducer Pair (22, 23) of Run 19, 

Phase I, cy = 14°, (3= 0°, M = 0.925, 6 /6 = 5°/12°, 6 = -10° 

(Concluded) ^ 


158 


FREQ. RANGE 
(Hz) 

RECORD PROCESSED 
( sec ) 

NUMBER OF 
RECORDS USED 

TOTAL PROCESSED 
DATA TIME (sec) 

2-100 

2.0 

15 

30 

100-800 

0.2 

16 

3.2 

800-4500 

0.1 

16 

1.6 


Note that the frequency ranges stated above were the test frequency scales. 
The power spectral density plots for all transducers (Nos. 14, 17, 22-25) cor- 
responding to Q/ = 0°, 8 °, 10°, 14° are presented in Figures 83 through 86 , 

respectively. All scales used in Figures 80 through 86 were modified to conform 
to the equivalent flight test condition of the F-5A, 

Referred to Figure 83, the high PSD levels of transducer Nos. 24, 25 at 
= 0 ° indicated possible shock-induced flow separation on the tail lower 
surface due to the negative tail deflection, 5 ^ = -10°. Good spatial correlation 
(Figure 58) with - 0.3156 at t' = .00955 sec (p^ = maximum correlation coeffic- 
ient, t' = time delay) between the two transducers was observed. However, the 
different auto— correlation function shapes indicated zero or minimum convection 
between transducer pair (24,25). 

As angle-of-attack increased to 8 degrees (Figure 84), moderate higher PSD 
levels on the upper tail surface (as well as on the upper wing surface) indicated 
the wake flow effect. On the other hand, the slightly lower PSD levels of trans- 
ducer Nos. 24, 25 (lower tail surface) implied that the severity of tail buffet 
was reduced by increasing the angle— of— at tack. 

At = 10°, the tail section chord was parallel to the mean flow direction. 

At this time, the dynamic pressures on the tail were caused essentially by the 
wake flow, and the rms coefficients for all locations on the tail surface were 
in the order of 0.1. Compared to = 8 ° case, the PSD levels showed moderate 
increase for all transducers. This increase was more pronounced for the fre- 
quency range up to 100 Hz. The approximately 15 dB increase in the PSD level 
of transducer No. 14 (Figure 85 ) signaled the local flow separation (also see 
Figure 79). Good spatial correlation existed between transducer pair (14, 22) as 
shown in Figure 81 indicating strong wake flow condition. The maximum value of 
the major peak of the correlation coefficient was 0.296 at t' = No clear 
convection mechanism was observed. Poor correlation existed between transducer 
Nos. 14 and 24. However, good correlation with p = 0.497 at < 7 ’ = 0.0015 sec. 
existed between transducer Nos. 24 and 25. 
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Beyond q. = 10°, the rms pressures on the lower surface of the tail de- 
creased, while those on the upper surface increased. At (<y = 14°, both trans- 
ducer Nos. 14 and 17 were under the separated flow. Spatial correlation 
between transducer pair (14,22) was weak, while transducer pair (22, 23) were wel 
correlated with = 0.5526 at t'= 0.0023 sec (Figure 82). Both transducers 
had similar PSD functions (Figure 86 ). The lag time suggested a convection 
speed of 236m/sec in the outboard direction, (see Reference 1, p. 109). Trans- 
ducer pair (24, 25) were also well correlated with = 0.4197 at = -0,0015 sec, 
but substantial difference in PSD pattern was observed. 

Additional processing of the wake flow and tail surface pressure data 
was performed with special attention directed toward the establishment of 
an overall picture of the wake flow pattern, as well as the effect of the 
tail surface setting under otherwise identical flow conditions. In Figures 
87 through 90 > the dynamic pressure coefficients measured at transducer Nos. 

14, 18, 22-25 of four cases of Phase I data corresponding to variable tail 
surface settings are plotted vs. angle-of-at tack . Transducer Nos. 14 and 
18 were located on the trailing edge flap, upper surface of the right wing 
at 47 percent and 33 percent semi-span locations, respectively. Transducer 
Nos. 22 and 23 were located on the upper surface of the horizontal tail 
and transducer Nos. 24 and 25 on the lower surface. Transducer Nos, 22, 

24 were located downstream of No. 14 while transducer Nos. 23,25 were located 
downstream of No. 18. The test configurations and the run conditions 
of Figures 87 through 90 are described below: 


FIGURE NO. 

TEST PHASE 
RUN NO. 

CONFIGURATION 

NO. 



M 

87 

1-8 

2 



50/120 

0 ° 

0.925 


1-19 

3 



5°/12° 

- 10 ° 

0.925 

88 

1-39 

9 



0 °/ 0 ° 

0 ° 

0.925 


1-49 

11 



0 °/ 0 ° 

- 10 ° 

0.925 

89 

I-I3 

2 



5°/12° 

0 ° 

0.750 


1-22 

3 



5°/12° 

- 10 ° 

0.750 

90 

1-43 

9 



0 °/ 0 ° 

0 ° 

0.750 


1-52 

11 



0 °/ 0 ° 

- 10 ° 

0.750 

Referring to 

Figures 87 through 90, the 

development 

of the wake flow pattern. 

the effect of 

the tail surface settings 

, and 

the horizontal 

tail buffet phena 

menon could be detected by the trend of 

the dynamic 

pressure coefficients 
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r 


as a function of angle-of-attack. In plots (a) of these figures 

(transducer Nos. 14 and 18), asterisks are used to indicate the angle-of- 

attack setting a * above which the C— values increased in such a manner as to 
s p 

indicate local flow separation. It may be seen that 


~ 5° 

for 

transducer 

No. 

18 

(337. 

semi-span) 

~ 7° 

for 

transducer 

No. 

14 

(477. 

semi-span) 


Furthermore, the solid lines denote the 6, = 0° cases and the dotted 

h 

lines denote the 6^ = -10 cases. In plots (b) (transducer Nos. 22, and 23, 
tail upper surface), asterisks are used to indicate the initial angles-of-attack 
Ck'^ ( ^ c^g) beyond which the horizontal tail surface flow pattern was influenced 
by the wake flow during wing buffet, where the rms dynamic pressures increased 
substantially. Wake flow effect on the horizontal tail upper surface could be 
traced for both 5^^ - 0° and -10° cases. For the configurations with 0 degree 
horizontal tail surface setting, the very high rms dynamic pressures of trans- 
ducer No. 22 at high angles-of-attack indicated the horizontal tail top surface 
buffet. The extremely high C- value (0.298) of Run 8 at O' = 16° (Figure 82b) 
was one of the highest rms dynamic pressures recorded in the test program. 

The flow patterns developed on the horizontal tail lower surface as 
shown in plots (c) of Figures 87 through 90 were quite different from those 
on the upper surface. For the configurations with zero degree of horizontal 
tail setting, high rms dynamic pressures were observed at low angles-of-attack 
(o' < 4° ) when the leading edge and trailing edge flaps of the model were 
deflected (solid lines in Figures 87c and 89c). It can be visualized that the 
high dynamic pressures on the tail lower surface were caused by the wake 
flow due to the deflected trailing edge flap. Otherwise, no wake flow effect 
due to wing buffet (except for transducer No. 24 of Run 13 in Figure 89c) 
was observed. 

For the configurations with -10 degrees horizontal tail setting (8^), the 
high C- values of transducer Nos. 24 and 25 at lower angle-of-attack settings 
(dotted lines in plots c of Figures 87 through 90 ) indicated shock-induced 
flow separation on the tail lower surface due to the negative tail deflection. 

As the angle-of-attack increased, the combined effects of tail buffet and wake 
flow due to wing buffet made the rms dynamic pressure levels of transducer 
Nos. 24 and 25 even higher for the configurations with deflected leading 
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edge and trailing edge flaps (Figures 87c and 89c). However, C— values 
dropped when the angle~of-attack setting exceeded 10 degrees where the 
negative horizontal tail setting cancelled the angle-of-attack setting. 

For the configurations with retracted leading edge and trailing edge flaps 
(Figures 84c and 86c) , the effect of wake flow due to wing buffet to the 

tail lower surface was minimum. The reduction in C- for transducer Nos. 

P 

24, 25 at very high (y (>12°) and 6^ = -10° also may be contributed to 
the wing downwash effect which tended to steer the wake flow downward 
and away from the horizontal tail. 

4.8 Correlation Study of Pressure Data Under Sideslip Conditions 

In order to further examine the sideslip angle effect on the buffet pres- 
sure distributions, the dynamic pressure data of Run 11 and Run 12, Phase II, 
were processed. Both runs featured the same model configuration except for the 
nominal sideslip angle which was +8 degrees for Run 11 and -8 degrees for Run 12, 
The configuration featured wing-tip-mounted Sidewinder missiles. The leading 
edge and trailing edge flap angles were 5°/12°. The aileron and horizontal 
tail surface settings were 0° and -10°, respectively. The Mach number was 
0.925. 

The majority of the pressure transducers were located on the right-hand 
side of the wing and the horizontal tail. Taking advantage of the model 

symmetry along the vertical plane, a comprehensive pressure distribution may 

be synthesized by making use of the pressure data of two sideslip angles (+8°). 

For instance, for Run 11, g = 8°, the leeward (left-hand) side pressure data 

were substituted by the leeward (right-hand) side pressure data of Run 12, 

3 = -8° . 

The rms pressure coefficient (C_) data for the two runs mentioned above 

P 

were processed. The initial development of the flow separation for Run 11 
was similar to those for Run 6, Phase II (Figure 34) which featured the 
same model configuration except that ~ However, up to q/ = 4° no flow 

separation was detected on the right wing upper surface. At = 6° shock- 
induced separation appeared at 617o and 477o semi— span, 607o chordwise positions 
(transducer Nos, 11 and 13). At o' ~ the separated flow region was extended 
outboard toward the tip. At = 10°, a major expansion of the separation bubble 


170 



took place. Except for the extended leading edge flap and the wing root area, 
most of the top surface of the right wing was in the separated flow region 
(Figure 91). As ^ reached 12 °, the right wing upper surface became completely 
separated. On the other hand, the high rms dynamic pressures on the horizontal 
tail lower surface at lower angle-of-attack settings ( 0 - 8 °) indicated the shock- 
induced flow separation due to the negative tail deflection = - 10 °. 

For Run 12, 0 = - 8 ° case, the flow patterns developed on the right wing 
upper surface (equivalent to the left wing upper surface for Run 11, 0 = 8 °) 
were quite different from those corresponding to Run 11. At = 4° local flow 
separated at 337» semi-span, 827> chordwise position (transducer No. 18) and also 
at 857, semi-span, 457 chordwide location and its downstream region (transducer 
Nos. 3, 4, 5). At Q- = 6 ° approximately 707, of the wing upper surfaces were in 
the separated flow region. At = 8 ° the flow was separated except for. an area 
on the extended leading edge flap. For a ~ 10°, the flow on the entire wing 
upper surface became separated (Figure 91). The rms dynamic pressures on the 
horizontal tail lower surface at low angle-of-at tack settings also indicated 
shock— induced flow separation. However, the C— values on the leeward side of 
the horizontal tail was 307, lower than that for the windward side at = 0°. 

In observing the development of the flow patterns on both the windward and lee- 
ward sides of the aircraft under sideslip, the following conclusions may be 
drawn: 

On the windward side, flow separation originated at the center portion 
of the wing semi- span, then extended outboard toward the tip, and then extended 
inboard toward the root section. On the leeward side, flow separation originated 
near the tip and extended toward the center. 

The flow on the leeward side wing upper surface became separated with 
lower angle-of-at tack settings as compared to the windward side. In general, 
it took approximately an additional 2 ° to develop an equivalent separated flow 
pattern on the windward side of the wing upper surface in comparison with that 
on the leeward side. Both the wing leading edge sweep-back and the fuselage 
wake effect are considered to be the leading contributing factors to the differ- 
ence in flow separation development on both sides of the wing with sideslip. 
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On the test model, three pressure transducers (Nos. 26 to 28) were installed 

on the left wing top surface at location symmetrical to transducers No. 2, 4, 

11 on the right wing top surface. In case the model featured perfect aerodynamic 

symmetry about its vertical plane, then the pressure data of transducer Nos. 

26-28 under a sideslip condition should be identical to those of transducer Nos. 

2, 4, 11 respectively when the sideslip angle was changed from positive to 

negative, or vice versa. To check the validity of this assumption, the C- 

P 

coefficients based on transducer Nos. 26 to 28 for 3=8° case were overplotted on 

the right wing C- data for g = -8° case in small circles. Figure 91 (a). Similarly, 

the C— coefficients based on transducer Nos. 26-28 for o = -8° case, were over— 

P ” 

plotted on the C- data of the right wing for g = 8° case. Figure 91(b). Allowing 

for minor irregularities in model configuration and variations in test conditions, 

the closeness of the left and right wing rms pressure data for equal and opposite 

sideslip angle (+8°) as indicated in Figure 86, seemed to justify the assumption 

of model symmetry and thus the feasibility of transferring the right wing pressure 

data to the left wing with the sign of the sideslip angle reversed. 

In order to study the pressure spatial correlations under the sideslip 
condition, the measured dynamic pressure data corresponding to = 10° were 
digitized. The pressure data digitized and processed were transducer Nos. 1, 3, 

4, 12, 14, 17, 22, 24, 27, Their locations are given below. 


T ransducer 

No. 


Surface 7 

o Semi-Span 7 

o Chord 

1 


RH 

Wing Upper 

85 

10 

3 


RH 

Wing Upper 

85 

40 

4 


RH 

Wing Upper 

85 

60 

12 


RH 

Wing Upper 

47 

44 

14 


RH 

Wing Upper 

47 

80 

17 


RH 

Wing Upper 

33 

42 

22 


RH 

Tail Upper 

47 

40 

24 


RH 

Tail Lower 

47 

40 

27 


LH 

Wing Upper 

85 

60 

Transducer 

Nos . 

22, 24 were 

located downstream of 

No. 14, and No. 

24 was 


directly under No. 22. Transducer No. 27 was located on left wing with the 
same relative position as No. 4 on right wing. The pressure data were digitized 
with a rate of 15,625 samples/sec and the corresponding folding frequency was 
7,812 Hz. A low pass analog filter with a cutoff frequency of 7000 Hz was used 
during the digitizing process. Applying the digitized data, spatial correlation 
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between transducer pairs were performed using a spectral analysis computer 
program. A time span of 0.131 second was used. The maxium time lag used in the 
cross -correlation computation was one-fourth the time span. The corresponding 
frequency increment was 15 Hz (approximately 2 Hz increment on the equivalent 
aircraft scale). 

The cross-correlation coefficient and the coherency function plots gener- 
ated by the computer program showed poor correlation between all transducer 
pairs except (1, 3) and (22, 24) for both p = 8 ° and p = — 8 ° cases. Typical 
plots are presented in Figures 92-98. The frequency scale was modified to 
conform to the equivalent flight test condition of the F-5A. In order to show 
the general trend of the coherency function in a linear frequency scale, the 
frequency range shown in Figures 92-98 was limited up to 157 Hz on the equivalent 
aircraft scale (approximately 1,200 Hz on the wind tunnel test frequency scale). 
The definition of cross— correlation coefficient and coherency function, may be 
found in Reference 1. 

Figure 92 showed the coherency function and the cross-correlation co- 
efficient distribution for transducer pair (1, 3) for 0 = 8 ° case. Under this 
test condition, the flow at transducer No. 1 was unseparated while No. 3 was 
approximately under the shock interaction line. The figure showed moderate 
correlation between the pair for frequency up to 55 Hz. The maximum cross- 
correlation coefficient was 0,36, For 6 = - 8 ° case. Figure 93, both trans- 
ducers were in the separated flow region. Pressures were rather well correlated 

with D = 0,7 and convection from No, 1 to No, 3 could be detected. The -spatial 
o 

correlation between transducer pairs (4, 27) and (22,24) are presented in 
Figures 94 - 97 . Figures 94,95 yielded the information about the degree of 
correlation between pressures on right and left wing upper surfaces. No 
significant correlation was detected. Figures 96, 97., showed that the degree 
of correlation of the pressures on upper and lower surfaces of right horizontal 
tail (transducer Nos. 22, 24) was almost as strong as that of two adjacent wing 
pressure stations for both g = 8 ° and g = - 8 ° cases. The cross-correlation 
function plots had a high peak at t = 0 and low level values throughout the 
remaining time scale. The negative peak value was confirmed by the phase 
difference of 180° for the two stations as shown in Figure 98. Furthermore, 
the PSD shapes for transducer Nos. 22, 24 were similar. All the observations 
gave clues as to the wake flow makeup near the tail surface as the tail 
section chord was parallel to the mean flow direction. The average phase 



z o iH r* rn DO o r> I i/» wi o X <*> zot-t-tozc"^ -<r>znxmio^ 


1-0 1.15 


F-5A WIND TUNNEL BUFFET TEST RUN 2,11 ID 79 ALP-10 


1.000 


0.800 


0.600 


O.MOO 


0.200 


.200 


.0 


Figure 


PRESSURE NO. 3 / PRESSURE NO. 1 




92. Spatial Correlation Plots for Transducer Pair (1, 3) of 
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Figure 93. Spatial Correlation Plots for Transducer Pair (1, 3) of Run 12, 
Phase II, 0 ^= 10°, P= -8°, M = 0,925, 6^/6^ = 5°/12°, 6^ = -10" 
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Phase II, a = 10°, P = -8°, M = 0.925, 6 /be = 5°/12°, 6, = -10° 
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Figure 96, Spatial Correlation Plots for Transducer Pair (22, 24) of Run 
11, Phase II, ^ = 10*^, P= 8° » M = 0.925, 6 /5^ ~ 5 /I2 , 6^ 
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Figure 97, Spatial Correlation Plots for Transducer Pair (22, 24) of Run 12 
Phase II, O' = 10°, P = -8°, M = 0.925, 6^/6^ = 5°/12°, 6^ = -10° 
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Figure 98, Modulus and Phase Angles of Pressure Cross- Spectrum for Trans- 
ducer Pair (22, 24) of Run ll, Phase II, O' = 10°, P = 8 , 

M = 0.925, 6 /6. = 5V12°, 6^ = -10° 
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angle difference of 180° seemed to be generated by the induced pressure due to 
random tail motion. 

4.9 Non-dimensional Power Spectra of the Dynamic Pressures 

In phase II, Runs 10 and 14 were performed under almost identical conditions 
Specifically, the test conditions are as follows: 

Model Configuration 3, - 5°/12° 

5 = 0°/0°, 5 = _ioo 

a h 

B = 0° 

M = 0.929 

= 7.447 X 10^/m (2.27 x 10^/ft) 

P = 29.90 KN/m^ (624.5 psf) 

P^ = 52.19 KN/m^ (1090 psf) 

Q = 18.05 KN/m^ (377.0 psf) 

For cy = 8°, a number of dynamic pressure data were converted into power 
spectra. The spectra were plotted in a non-dimensional form so as to compare 
with non-dimensional power spectra data collected by Dennis G. Mabey (Reference 
4) according to the convention orginated by T. B. Owen (Reference 7), In both 
references, the data were plotted as: 

/nF(n) vs n = fL/V 
where f = frequency, Hz 

L = a characteristic length (the bubble length in Reference 4) 

V = the freestream velocity 

2 2 

F(n) = contribution to < P > /q in frequency band Af 

P = the fluctuating pressure 

Q = q = the freestream dynamic pressure 

In Figures 99 to 101, the non-dimensional pressure spectra for transducer 

Nos. 18, 22, 23 are plotted for F-5A model corresponding to = 8°, Run 14 » 

Phase II described previously. The coordinators are: 

2 

0f/Q vs k = o)b/U 

where 0 = the power spectral density 

00 = circular frequency = 2Trf 

b = a characteristic length (the mean aerodynamic chord) 

U = V = the freestream velocity 
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PRESSURE NO. 23 



k = u)b/u 


Figure 101 Non-Dimensional Pressure PSD of Transducer No. 23 Recorded at 

a = 8°, Run 14, Phase II, p = 0°, M = 0.929, 6 / 6^ = 5° f 12\ 6^ = -10° 
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Compared with the British convention, it is seen that: 

0f/Q^ = nF(n) 

k = 2rt(b/L)n 

SO that the spectral plots of Figures 99 to 101 can be compared readily with the 
British data shown in Reference 4. In the above, the ratio (b/L) is the ratio 
between the aircraft mean aerodynamic chord and the separation bubble length. 
Except for the PSD peak at k = .05 for transducer No. 18 corresponding to the 
first symmetrical wing bending mode, the general patterns of the non-dimensional 
PSD' s of the F-5A data are qualitatively comparable to those compiled in 
Reference 4. 



Section 5 


ANALYTICAL PROGRAM - SCALE MODEL STABILITY PREDICTION 
ALLOWING FOR AEROELASTIC EFFECT 

The scale model used in the wind tunnel test is a structure with a 
flexible support. The spring constant associated with each component of 
the flexible support may be determined as the product of the square of the 
measured natural frequency (in rad/sec) and the total mass or mass moment- 
of~inertia of the scale model associated with that degree of freedom. For 
flexibly-supported structures, the equations of motion in the frequency domain 
ay be written as: 


^11 ^12 
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^12 



'ho 

. + i(« 



L 21 22 S 
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L 21 

22 J 



where K, C, M are the stiffness, damping, and mass matrices of the structure; 
h and f represent the Fourier transforms of the deflections and external forces, 
respectively. In Equation ( 8 ), subscript 1 refers to the control points and 
subscript 2 refers to the reference point where the flexible support is attached. 
K is a diagonal matrix whose elements represent the spring constants associ— 

O 

ated with various degrees of freedom of the flexible support. It is noted that 
an infinite spring constant implies a rigid constraint, while a zero spring 
constant implies a free flight condition. The values of K in Equation ( 8 ) are 
assumed to be finite, and that: 


K . 

SI 



si 


(9) 


In Equation ( 9 ), M is the total mass or mass moment— of— inertia of the structure 

si 

in the direction of the i-th degree of freedom, and m . is the measured nautral 

SI 

frequency in rad/sec. To simplify the calculation procedure, the mass of the 
structure is assumed to be concentrated at the control points. 


The governing equation for an undamped natural vibration of flexibly- 
supported structure is simple: 


*^11 *^12 
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where x is the mode shape and u)j, is the resonance frequency. Equation (10) may 
be simplified by introducing a control point rigid transformation matrix T whose 
columns are composed of rigid body modes associated with the unit movement of 
each degree of freedom of the flexible support, respectively. The matrix T 
satisfies the following relationship; 


■'^11 '^ 12 ‘ 


n 

_^21 *^ 22 . 


LiJ 


= [ 0 ] 


( 11 ) 


K^2 = T 


K22=T 


(12) 


Since, for K = 0, the structure is free to move as a rigid body in the flexibly- 
supported directions, substituting Equations (9) and (12) into Equation (10) gives: 




K^lT + Mg ^ 





(13) 


and the eigen equation is then: 


- 

FT 2 

-1 .t,T 

K - K T 
11 11 

T + 

T 



(14) 


Let X be a modal matrix whose columns are composed of modal shapes x. The 
orthogonality conditions between modes lead to the following diagonal matrix 
expressions for modal masses and modal stiffnesses; 




(15) 


With the available modal data, the structure deflections may be represented as: 


hi 

(■’2 

where is the modal amplitude matrix. Premultiplying Equation (1) by X^ 

and using Equations (15) and (16) yields the modal equation of motion in the 
following form; 
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provided that the damping matrix Cc] satisfies the following conditions; 
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(18) 


where g„ is the modal damping coefficient matrix. 

F 

The instability characteristics along the roll axis of the flexibly-supported 
model may be determined by the use of the dynamic equations applied in stab- 
ility and control studies. The lateral stability equations of a rigid air- 
craft in terms of the stability derivatives derived with respect to a body 
axes system with the origin at the center of gravity and x-z plane the plane 
of symmetry, may be written as follows : 




2 ir 2 Zp 


mU (B + cos a r - sin a p) = (C UB+^C r+^C p) 
o o U yB 2 yr 2 yp ^ 


; = asb 


I r - I p ,, 
zz xz U 


(C U B + T C r + ^ C p) 
n3 2 nr 2 np ^ 


(19) 

( 20 ) 
( 21 ) 


where p, B » r are the roll rate, the side-slip angle, and the yaw rate, 
respectively, m is the mass, 1^^> moments of inertia about the 

roll and yaw axes, respectively; is the product of inertia. U is the flight 

speed; is the initial angle-of-attack; b is the span of the aircraft and S is 
the wing area, P is the air density and q is the dynamic pressure. The other 
expressions in Equations (19-21) are the non-dimensional aerodynamic derivatives 
that are the functions of aircraft configuration, the altitude, the angle-of- 
attack, the Mach number, and the Reynolds number. For a given flight condition, 
the solution of the eigen-equations (19-21) yields the aircraft modal data, includ- 
ing the Dutch roll oscillation which is of special interest in aircraft buffet 
studies. In applying the dynamic Equations (19-21) to the flexibly-supported 
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scale^ model, proper spring and damping terms representing the supporting systeu 
reactions as well as a number of fundamental wing bending-torsion modes were 
added. The dynamic equations in the matrix form are described below. 
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where a are the modal amplitudes. In Equation (22) ^ subscrip t r refers to the 
supporting system modes and subscript f refers to the structural flexible 
modes. The expressions of the modal spring, damping and mass terms (K, C, 
M)^ ^ are given in Equations (15)and (18). A^ ^ represents the generalized 

aerodynamic matrices. The eigen-value A is in general a complex number 


A = iw (1 + iC) 


(23) 


where (jo is the angular frequency and ^ the damping ratio. The criterion 


^ > 0 


indicates that the motion is stable. 

Let 

{a^} = ^ roll angle, side deflection, yaw a ngle ^ (24) 

then the matrix A^^ in Equation (22) may be approximated in terms of the 
nondimens ional aerodynamic derivatives by means of Equations (19 to 21) 
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In predicting the instability characteristic of the flexibly-supported model, 
the scale model modal data were .generated under the following conditions: 

1. The weight and the masQ moment of inertia about roll axis of the model 
match their measured values. 

2. The natural frequencies, corresponding to the spring support modes and the 
first wing bending (- torsion) mode match their measured values. 

The F-5A basic aerodynamic stability and control data (Reference 12) were 
used to determine the aerodynamic derivatives, which are functions of altitude, 
Mach number and angle-of-attack, under various test conditions. 


Using the above formulation, the analytical computation of the instability 
characteristics of the flexibly-supported scale model for a = 0°, 4°, 8°, 

12® was performed for the following configurations: 


Mach Number 

Free stream speed 

Free stream dynamic pressure 

Equivalent altitude 


0.925 
299 m/sec 
18.37 KN/m^ 
10,668 m 


0.75 

249 m/sec 
15.9 KN/m^ 
7,772 m 


Three supporting system modes and two structural flexible modes were used in 
the analysis. The damping coefficient corresponding to the roll motion was 
0.024, the lowest damping level used during the tests. 


The calculated instability characteristics along the roll axis of the 
model are presented in Figure 102. Since the natural frequencies generated 
by the eigen -value problem were so close to their undamped values, only the 
damping ratio C (Equation 23) is plotted vs. angle-of-attack. 


Referring to Figure 102, the damping for the wing rock type oscillations 
of the scale model decreases rather rapidly with increasing angle-of-attack 
beyond 8 degrees. It can be visualized that, for each configuration, there 
was a critical angle-of-attack a* at which the model started to oscillate in 
a random manner. For the scale model tested, wing rock appeared at approx- 
imately a* = 8° to 12® for M = 0.925, depending on the test configurations. 

For M = 0,75, the values of qj* are higher in comparison with those correspond- 
ing to M = 0,925 cases. This trend is consistent with the analytical predic- 
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tion as shown in Figure 102. In a free flight condition, near zero or no mechani- 
cal damping exists. For this case, the zero damping level will be moved to the 
C — 0,012 line. As a result, the - q/ curve of Figure 102 indicates the 
aircraft becomes unstable at a =- 14® for M = 0,925. 

A similar approach could be applied to predict the aircraft instability 
characteristics in free flight conditions by dropping the K^., terms in 
Equation (22). However, it must be remembered that the results of the wind 
tunnel test can only be applied to free flight stability considerations in a 
qualitative manner. This is because the yaw and sideslip responses, which are 
important in Dutch roll, for free flight were not adequately modelled in the 
wind tunnel test. In the following, the results of predicted Dutch roll motion 
of F-5A are presented. 



M = 0.925, h 

= 10,668 m 

M = 0,75, h = 

7,772 m 


f (Hz) 

C 

f (Hz) 

C 

a 

= 0° 

0.412 

0.1009 

0.34 

0.1087 

a 

= 4° 

0.456 

0.1161 

0.372 

0.131 

a 

= 8® 

0.536 

0.1339 

0.44 

0.1476 

a 

= 12® 

0.619 

0.0805 

0.493 

0.1115 


Based on the above data, the following conclusions are drawn: 

1, The highest aerodynamic damping occurs at a ^ 8°. Beyond this point, 
damping decreases rather rapidly indicating wing rock instability at 
higher oi angles. 

2, In free flight condition, the predicted wing rock occurs at lower angle-of- 
attack when M = 0,925 as compared to the M = 0,75 case. This conclusion 
is in line with the test results and the analytical prediction for the scale 
model. 
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Section 6 


CONCLUSIONS 

In the report, the wind-tunnel test program, the processed data, and the 
correlation results between the wind tunnel data and the buffet flight data 
are presented. Also included are the results obtained based on the parame- 
tric studies of the wind-tunnel data (by varying either the model configura- 
tion or the test condition). In the following, the key conclusions and ob- 
servations are presented. 

1. The separation region was developed over the upper wing surface with in- 
creasing angle— of— attack. The development of the separated flow over a 
"substantial" part of the outer wing occurred within a relatively narrow 
angle-of-attack range (A o' ~ 2°). For tests conducted at a fixed angle-of- 
attack in this range, the dynamic pressures were found to be highly non- 
stationary, particularly in the shock/ bound ary layer interaction region. 

2. In general, the development and expansion of the flow separation region on 
the model wing surface were similar to the results acquired during a previous 
buffet flight test program of the F-5A aircraft. At M = 0.925, the shock 
boundry was initiated near the trailing edge of the wing-tip region and grad- 
ually expanded toward the inboard and upstream region with increasing angle- 
of-attack. For M = 0.75 the flow separation was leading edge induced, 

3. In comparison with the zero sideslip angle case, the effect of positive 
sideslip angle (for the windward side of the wing semi-span) tended to move 
the shock and shock-induced flow separation downstream and to create a favor- 
able flow condition in its wing-tip region which retarded the local flow 
separation. On the leeward side of the wing semi-span (negative sideslip 
angle condition), flow separation seemed to originate at the wing-tip region. 
Furthermore, the flow on the leeward side wing upper surface became separated 
with lower angle-of-attack settings (approximately 2®) as compared to the 
windward side. 
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4. As compared to the cases with tip-mounted missiles, the lack of tip- 
mounted missiles and their trailing vortex flow seemed to aggravate the 
situation and make the shock location and the resulting flow pattern less 
stable. However, the tip— mounted missiles had an adverse effect in pro- 
moting the flow separation at a certain critical angle. This effect was 
greatly diminished after the critical angle-of— attack was surpassed. 

5. The deflection of the horizontal tail (to -10®) tended to delay or alter- 
nate the shock-induced flow separation in that part of the wing direction 
upstream of the tail. Correspondingly, for a moderate angle-of-attack, the 
shock strength and the level of dynamic pressure behind the shock tended to 
be stronger in the outer half span of the wing as compared to the zero 
horizontal tail surface setting case. The effect of the tail surface setting 
became less prominent as the angle-of-attack was further increased, 

6. At moderate angle-of-attack, the deflected leading edge flap tended to push 
the flow separation boundary further downstream as compared to the case of the 
retracted flap. The extended trailing edge flap tended to keep the shock-induced 
flow separation region further upstream in the inboard region at a moderate 
angle-of-attack,. This can be attributed to the increase camber in the inboard 
portion of the wing semi— span for the extended flap cases. 

7. For angles-of-attack under the critical angle (the angle-of-attack at 

which the flexibly-supported model started to oscillate along the roll axis), 
the trends of the shock development and the expansion of the flow separation 
region were the same for the fixed and flexibly— supported models. However, 
substantial deviation in local pressures did exist. At the critical angle-of- 
attack (y'‘, the flexibly-supported model motion tended to make the dynamic 
pressure distribution more even, spacewise , as compared to the fixed model. For 
Q( > the difference in the dynamic pressures of the models tended to diminish 

as the angle-of-attack was further increased. 
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8. When maximum wing-rock motion occurred, the dynamic pressure peak levels 
corresponding to the natural roll frequency were extremely high for most 

of the transducers of the flexibly-supported model, in comparison with the 
corresponding pressure data for the regularly-supported model. 

9. For the configurations with 0 degree horizontal tail surface setting, 
the very high dynamic pressures on the horizontal tail upper surface at high 
angles-of-attack indicated the horizontal tail upper surface buffet. On the 
lower surface, high dynamic pressures caused by the wake flow due to the 
deflected trailing edge flap were observed at low angles-of-attack. 

10. For the configurations with -10 degrees horizontal tail setting, the high 
dynamic pressures on the horizontal tail lower surface at lower angles-of- 
attack indicated shock- induced flow separation due to the negative tail de- 
flection. As the angle-of-attack increased, the combined effects of tail 
buffet and wake flow due to wing buffet made the d 3 mamic pressure levels on 
the horizontal tail lower surface even higher for the configurations with 
deflected leading edge and trailing edge flaps. However, dynamic pressure 
levels dropped when the angle-of-attack setting exceeded 10 degrees where the 
negative horizontal tail setting cancelled the angle-of-attack setting. The 
effect of wake flow due to wing buffet on the dynamic pressures on the tail 
lower surface was minimum. For the configurations with retracted leading 
edge and trailing edge flaps, the reduction in the dynamic pressure levels 

on the horizontal tail lower surface at very high angles-of-attack also may 
be contributed to the wing downwash effect which tended to steer the wake 
flow downward and away from the horizontal tail. 

11. The trend of development and expansion of the flow separation region were 
similar for the wind-tunnel model test and the corresponding flight test. In 
general, measured in angle-of-attack, it took the scale model approximately an 
additional 2° to 4° to develop a similar separation flow pattern on the. top wing 
surfaces. The Reynolds number effect was less significant. For M = 0.925 
case, at a given spanwise location, the shock tended to stay closer to the 
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trailing edge on the scale model as against the flight test results. Certain 
factors were noted which were believed to contribute to the scale model flow 
behavior; (1) the aeroelastic effect of the full-scale aircraft caused an 
increase in the local angle-of-attack in the outer span of the wing; (2) the 
transient effect in a transonic maneuver of the aircraft where the angle- 
of-attack changed continuously; (3) F-5A wing has a thin airfoil section, 
namely NACA 65A004.8; (4) transition strips were installed on this F-5A scale 
model which should trigger the shock induced flow separation in an upstream 
location on the wing chord. 

12. For given test condition, the instability characteristics along the roll 
axis of the flexibly-supported model was predicted vs. angle-of-attack using 
the dynamic equations employed in stability and control studies. A similar 
approach may be applied to predict the aircraft instability characteristics 
in a free flight. With additional refinements in the analytical formulation, 
the approach can be extended to include the steady state sideslip angle as 

a parameter in the wing rod. stability prediction. 

13, The flexibly-supported model is important in simulating the wing rock 
oscillations in a transonic maneuver. Based on the pressure and response 
data acquired when the model was near stationary and when it was excited by 
flow to oscillate in roll, the following conclusions may be drawn. 

A. Wing rock was generated by a limit cycle mechanism due to the fluct- 
uating pressure changes on the wing top surface, especially near the 
wing tip region. In the cases observed, the major reason for the 
pressure changes was the motion induced Aa which alternately caused 
leading edge stall and recovery. The amount of energy input to the 
limit cycle, to be consumed by aerodynamic and/or mechanical damping, 
may be evaluated quantitatively. 

B. Based on kinematics, the roll oscillations acted together with 

6 to cause substantial Aa changes which contributed to the limit 
cycle. Wing rock amplitude was thus dependent on the steady state 
sideslip angle B. 
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C, With a sideslip angle 3> the upper wing surface on the windward side 
contributed the majority of the fluctuating pressures which caused 
wing rock. 

D, With AcLy the local fluctuating pressure coefficient followed the ramp 
curve during leading edge stall and the consequent loss of lift. 
Dynamic overshoot was observed in both the lift loss and recovery 
phases . 

E, Based on the observations of the wind tunnel test data, a preliminary 
mathematical model of aircraft wing rock was developed which illus- 
trated the motion-coupled limit cycle mechanism. Since the limit 
cycle defining the energy input to the aircraft was dependent on a 
number of flight condition parameters, no uniqueness condition c^n 

be established at this time to define and predict wing rock. Never- 
theless, it is known that wing rock was directly related to the lift 
loss (ramp) of the vs. a curve covering the flue tua ting-pressure- 
active region of the wing. Additional acquisition and processing 
of controlled test data as those described above may further refine 
and quantify the preliminary mathematical model of aircraft wing rock 
described in this report. 


14. In actual transonic maneuver flight tests of F-5A aircraft, two distinct 

points were identified as the buffet onset and wing rock with increasing c^. 
The corresponding model tests at fixed ly conditions did not reflect the 
buffet onset point because of the substantial mechanical damping of the 
model. The critical or'' angle, at which the model became unstable and 
oscillated in roll, was in most cases within ± 2® of the wing rock o’ 
recorded in actual flight tests. 
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